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Abstract 

We consider the semiclassical asymptotics of the sum of negative eigenvalues of the 
three-dimensional Pauli operator with an external potential and a self-generated mag- 
netic field B. We also add the field energy f3 J and we minimize over all magnetic 
fields. The parameter /3 effectively determines the strength of the field. We consider 
the weak field regime with /3/i^ > const > 0, where h is the semiclassical parameter. 
For smooth potentials we prove that the semiclassical asymptotics of the total energy 
is given by the non-magnetic Weyl term to leading order with an error bound that is 
smaller by a factor h^~^^ , i.e. the subleading term vanishes. However, for potentials with 
a Coulomb singularity the subleading term does not vanish due to the non-semiclassical 
effect of the singularity. Combined with a multiscale technique, this refined estimate is 
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used in the companion paper |EFS3j to prove the second order Scott correction to the 
ground state energy of large atoms and molecules. 

AMS 2010 Subject Classification: 35P15, 81Q10, 81Q20 
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1 Introduction 

An important problem in semiclassical spectral analysis is to determine the sum of the negative 
eigenvalues of a Schrodinger operator —h'^A — V{x) on L^(R'^), i.e., 

Tr {-h'^A - V{x))^. 

We will use the convention that X- = = min{a;, 0} when x is either a real number or a 
self-adjoint operator. It is well known that under appropriate integrability conditions on V 
the leading order term is given by the Weyl asymptotics 

Tr (-/i^A - = (27r/i)~^ /" /" {p^ - V{x))^dxdp + o{h~'^), h ^ 0, (1.1) 

and for smooth potentials the error term can be improved to 0{h^'^~^'^) (under some non- 
criticality assumption) by using pseudo-differential calculus. In other words, the subleading 
term in the semiclassical expansion in powers of h vanishes. We remark that with more 
elementary methods and under less regularity assumptions on V, for a local version of this 
problem 

Tr [iPi-h^A - V{x))ij]^, e C^{R'^), 

the error term has been shown to be 0{h~'^^^^^^) in Theorem 12 of |SS] (we will recall it in 
Theorem 13. 3p . This bound was extended to the relativistic case in [SSSj . 

A generalization of this problem is to consider not only a potential V but also an exterior 
magnetic field B = V x A generated by a vector potential A. The corresponding magnetic 
Schrodinger operator is {—i\/ + A)'^ — V{x). A further generalization is to consider the particles 
as having spin— | and introduce the magnetic Pauli operator [cr ■ (— zV + A)]^ — V, where in 
d = 3 dimensions cr = (o"i, (72, (J3) denotes the vector of Pauli matrices. For simplicity we 
will consider the d = 3 dimensional case only and we denote both the Schrodinger operator 
{—ihW + AY and the Pauli operator [cr ■ [—ihV + A)]'^ by Th{A). Our analysis will be carried 
out in the more complicated case of the Pauli operator, analogous but easier results can be 
proved for the Schrodinger case as well. 
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Much work has gone into understanding the semiclassical asymptotics of the sum of neg- 
ative eigenvalues, i.e., the asymptotics for small /i > of 



Tr {Th{A) - V{x)).. 

It is well known that under appropriate conditions on A and V the leading behavior as h 
tends to zero is given by the Weyl asymptotics (11. ip and note that the limit behavior is non- 
magnetic, i.e. fixed magnetic fields do not influence the leading order semiclassics. If the 
magnetic field is rescaled, B fiB, and the coupling constant fi increases along with the 
/i — )■ semiclassical limit at least as /i > h~^, then magnetic fields become relevant even 
in the leading term. Most work in this direction has been carried out with a homogeneous 
magnetic field |Sobt ILSYlt ILSY2] with some generalization to an inhomogeneous one |ESt H] 
but always subject to regularity conditions on the field. 

In this paper we will address a related and equally important issue, namely the case when 
the magnetic field is not a fixed external field, but the self-generated classical magnetic field 
generated by the particles themselves. The vector potential A will be optimized to minimize 
the total energy consisting of the sum of negative eigenvalues (corresponding to the ground 
state energy of non-interacting fermions) and the field energy 

j B^ = j iVxAp 

(we use the convention that unspecified integrals are always on w.r.t. the Lebesgue mea- 
sure). The problem we consider is thus to determine the energy 

E{P,h,V) = M 

for /3, /i > 0, where the infimum runs over all vector fields A G if^(]R^ M^); in fact mmimizmg 
only for all A G C^(M^;R^) gives the same infimum. Alternatively, in addition to A G 
if^(M^;R^), one could impose a gauge fixing condition, e.g., V ■ A = 0, see more details in 
Appendix A of [EFST] . 

Here /3 is an additional parameter setting the strength of the coupling of the particles to 
the field. Formally /3 = oo corresponds to the non-magnetic case; smaller /3 means that a 
larger optimizing magnetic field is expected. In a given physical system the values of h and fi 
are given, but as is standard in semiclassical analysis we leave them as free parameters. 

The Euler-Lagrange equation corresponding to the variational problem (II. 2p above is 

/3Vx5 = Ja, (1.3) 
where J a is the current of the Fermi gas, which in the Schrodinger case is 
Ja{x) = -Re [{-ihW + A)l(_<^,o](T/.(A) - V)] (x,x) 



Ti{n{A)-V)^ + P / |Vx A 



;i.2) 
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and in the Pauli case is 



Ja{x) = -Re [Tr ■ (-zW + A))l^_^,o]{n{A) - V))] (x, x). 

In other words the Euler-Lagrange equations are the non-hnear coupled Maxwell-Schrodinger 
or Maxwell-Pauh equations. 

Semiclassical resuhs with magnetic fields mentioned above assume that the field is regular. 
However, if the magnetic field arises as self-generated and thus determined internally via a 
variational principle, the sufficient regularity is not a-priori given. 

The first semiclassical result for the local problem with a self-generated magnetic field 
was presented in |ES3j . where the leading order asymptotics was shown to be given by the 
non-magnetic Weyl term in the weak field regime, > const > 0. The error term was by 
a factor of order h^^"^ smaller than the leading term (see Theorem 1.3 of [ES] for the lower 
bound; the matching upper bound was not explicitly stated in |ESj but it clearly follows by 
choosing A = 0). 

The main result of this paper is a substantial improvement of the error term to a factor 
of order /i^"*"^. This result can also be interpreted as showing that the subleading term in the 
semiclassical expansion in powers of h vanishes. 

Theorem 1.1. Let V G C^(M'^). There exist a universal constant e > such that for any 
fixed kq > we have 



Hm h"^-' 

/1-5-0 , l3h'^>Ko 



E{/3,h,V) -2{2nh)-^ J j [p'^ - V (q)] _dpdq =0 (1.4) 



for the Pauli problem. The same result (without the pref actor 2) holds for the Schrddinger 
problem as well. 

In fact, we can replace the infimum over all A in the definition of E{(3, h, V) with a good 
apriori bound on the field energy: 

Corollary 1.2. There exist universal constants e > and Kq > such that if Ah is a sequence 
of vector potentials satisfying 

lim/i^-"/3 /"iV X AT = 0, (1.5) 



then for any fixed smooth, compactly supported potential V , we have 



TT{n{Ah)-V). + (3 j \VxAh\'-2i2nhy^ j j [p" -V{q)\_dpdq 







lim h'^-' 

h^O , /3/i2>fto 

(1.6) 

for the Pauli problem. The same result (without the pref actor 2) holds for the Schrddinger 
problem as well. 
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In |EFS1] we have also analyzed the semiclassical behavior of E{P, h, V) in other regimes 
of the parameter /3. The main motivation for a precise second order asymptotics in the specific 
regime Ph"^ > const > is that Theorem 11.11 is the main technical input for the proof of the 
Scott correction term of the ground state energy of large atoms and molecules in the limit 
when the nuclear charge Z tends to infinity. We refer to Section 2.3 of |EFS1] for the precise 
statement on the Scott correction with a self-generated magnetic field and for an explanation 
on its connection with semiclassical asymptotics. The actual proof of the Scott correction is 
given in a separate paper |EFS3] . 



2 Localized models 

The semiclassical asymptotics is essentially a local issue. If the potential V has sufficient 
decay there will be no semiclassical contribution from infinity. It is an interesting question 
whether the decay needed for stability will ensure that there is no semiclassical contribution 
from inifinity, but at this point we do not have a general proof of this. To separate the local 
semiclassics from the issue at infinity we introduce local versions of the energy. 

Let ip he a smooth function with supp'?/' C -B(l), where in general B{r) denotes the ball 
of radius r centered at the origin. We will always assume that HV^Hoo < 1 and ip is identically 
1 in a neighborhood of the origin. Denote by ipr{x) = ip{x/r). For any r > define 

EriP, h, V) := inf [Tr [MTh{A) - V)4j,] _ + /3 [ | V x ^l^j , (2.1) 



where we minimize over all A G or, equivalently, over all A G C^(M'^; M^). 

Alternatively, we may again restrict to vector potentials A G -ff^(M^;M^) with V ■ A = in 
(12. ip without changing the value of the infimum. Without loss of generality, we can always 
assume that V is supported in B{r). 

Here we have localized the particles but not the fields. We can also do both. Let us first 
note that if V ■ A = then 

|VxA|2= / |V®A|2, (2.2) 



where the integrand of the right contains all first derivatives, i.e. |V ® = Yl^j=i l^iAjl"^. 
This identity is easily seen using the Fourier transform. As the local version of the field energy 
we will use the localization of the integral on the right. For any R> r > we thus set 

Er,R{P,h,V):=M\TT [ip,{Tf,{A)-V)^r]_ + f3 [ \V0A\^], (2.3) 

^ L Jb(R) ^ 

where we minimize over all A G H^{B{R); M^)^ or, equivalently, over all A G C(j"(M'^; M^). The 
localized field energy is not fully gauge invariant, we therefore cannot restrict attention to 
divergence free vector potentials. We are however free to add a constant vector to A. Notice 
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that the integral of |V ® is taken on a larger ball B{R) than B{r) which contains the 
support of ipr- By virtue of gauge invariance and fl2.2p we have 

Er = Er oo '■= lim Er R. 

R^oo 

By a simple rescaling we have 

Er,R{P,h, V) = r~'E,,n/r{rP,h,Vr), 

where Vr{x) = r'^V{xr). It is thus enough to analyse the case r = 1. 

The definition fl2.ip is physically somewhat better motivated than f l2.3p since it contains the 
energy of the magnetic field in the whole space and thus gives the correct magnetic interaction 
between the particles. The form (12. 3 p is however more useful if we want to localize all parts 
of the energy. 

Another version of the localized energy would be 

Elj,{P,h,V):=ini\TTi;^^[mA)-V]_ + P [ iV^Apl. (2.4) 

^ L Jb{R) ^ 

This has the disadvantage of being more complicated to calculate as it requires knowledge of 
the operator Th{A) — V on the whole space, e.g. it is not enough to know V only on B{r) or 
A only on B{R). It has the advantage of not causing localization errors. Note that since il) is 
assumed to be equal to 1 in a neighborhood of we may identify 

E(/3, h, V) = hm E;^(/3, h, V) =: E'^^^{(3, h, V). 

r— >cxD 

Finally we could also have defined the localized energy by introducing Dirichlet boundary 
conditions on the boundary of B{r), i.e., 

E^j,{{3,h,V) -.= ^^1 [{Th{A) - V)B(r),D]_ + 13 [ |V®A|2]. (2.5) 
^ L Jb(r) ^ 

We note that the leading order local semiclassical result (Theorem 1.3 in jESj ) was stated for 
Dirichlet boundary conditions. 

In the main part of this work we will mainly use the global energy E in (II. 2p or the localized 
energy Er^R in (12. 3 p but in the next section we will compare the local versions Er, Er^R, and 
E'r R- We will not discuss Dirichlet boundary conditions. 

2.1 Comparison between the different localized energies 

We first compare the energies Er^pt and E^-. 



6 



Lemma 2.1. There exists a universal constant Cq such that for all < r < R/2 we have 

ErAf^, h, V) < EriP, h, V) < Er,R{{l + Co{r/R)^)/3, h, V). (2.6) 
This result holds for both the Pauli and Schrddinger operators. 

Proof. The first inequality in fl2.6p is trivial since Er = -E'r,oo and Er^R is clearly an 
increasing function of R. 

To prove the second inequality we start from an approximate minimizing A G for the 
energy Er^R on the right hand side of fl2.6p . By subtracting a constant vector from A we 
may assume that A has average on the sphere |a;| = r. We may, moreover, assume that 
A minimizes the Dirichlet integral |V ® v4(x)pda; given the boundary value of A on 

|x| = r. If not, we could improve the trial energy of A by replacing it on the set r < |x| < R 
with the vector field that agrees with A for |a;| < r and minimizes Jj-^^^^^r |V ® A{x)\'^dx. As 
a trial vector field for the Er we choose the field A' defined on all of M.^ that agrees with A 
for |x| < r and minimizes the integral |V ® A'{x)\'^dx. The vector fields A and A' in 
the region |x| > r can be expressed in terms of the common boundary value of A on |a;| = r. 
From the lemma below we see that there exists a universal constant Cq such that 

f |V® < (1 + Co(r/i?)3) /" |V®A(x)pda;. 

Jr<\x\ Jr<\x\<R 

Note that the constructed A' is not necessarily divergence free. Since A and A' agree on 
|x| < r we obviously have 

[ |V X A'(x)pdx < /" |V® A'(x)pdx < (l + Co(r/i?)3) /" \V ® A{x)\^dx 

Jr3 JR3 J\x\<R 

and 

Tr [iPr{Th{A) - V)^r\_ = Tr [MTh{A') - V)iJr]_. (2.7) 
The second inequality of (12. 6 p follows from the last two observations. □ 
We now give the simple estimate needed in the previous proof. 

Lemma 2.2. Assume < 2r < R and let S{r) = {\x\ = r}. Given g G L'^{S{r)) with 
average 0. Let fi G H^{B{R) \ B{r)) and f2 G H^{M.^ \ B{r)) satisfy the boundary conditions 
fi\s{r) = f2\s{r) = 9\s{r) o,nd minimize the respective Dirichlet integrals 



B{R)\B{r) jR3\B(r) 

Then 

\Vf2\'<{l + Co{r/Rf) I |V/ir, 



/ 



(r) JBiR)\Bir) 

for some universal constant Cq. 
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Proof. By rescaling we may assume that r = 1. We know that /i is harmonic for r < |x| < R 
and satisfies Neumann boundary conditions on the boundary \x\ = R and that /2 is harmonic 
on |x| > 1 and tends to at infinity. We can write 

/ |V/ir = -/ gdrfu [ \Vf2? = -f gdrf2. 

Jb{R)\B{1) Js{i) Je3\B(1) Jsil) 

We expand g in angular momentum eigenfunctions g = J2'^o'^m=-e'^^mYem- Since g has 
average 0, we have aoo = 0. We also expand the harmonic functions 

fli^) = X](«l£mk|^ + burn\x\~'^~^)Yira{x/\x\) 
Im 

and 



f2{x) ='^h2im\x\ ^ ^Yim{x/\x\). 

im 

The equality of the boundary conditions at |x| = 1 implies 
The Neumann condition at |x| = R implies 
(Note all coefficients vanish for £ = 0). Hence 

htm={l+r\i+l)R-^'-^)hun.. 

We thus obtain 



'S{ 

and 

oo e 



/ gd,.h = 4vrV V(l-i?-2^"i)(l + r^(£+l)i?"2^-i)(£+l)|6i,^| 

"^^(1) i=l m=-l 



> -{1-3R-') / gdrf2. 

'5(1) 



□ 
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Now we compare Er^R and E^j^. Note that the only difference between Er^R and E^.^^ is 
that in fl2.4p ipr is outside the negative part. The following lemma shows that Er^R{f3, h, V) 
and E'j. h, V) are essentially equivalent up to a term of high order in h (to appreciate the 
next result, recall that Er^R{(3, h, V) is typically of order h~^). 

Lemma 2.3. There exist universal constants Co,Ci,C2 > 0, such that 

E',^^{(3,h,V)<Er,R{/3,h,V) < E:.^^{{l + Co{r/Rf + 6)f3,h,V) (2.8) 

+Cih~\\\V\\l^/l + r~') + C,h'6-'p-'{\\V\\t + r-') 

for all (5 > 0. The same result holds for the Schrddinger case, with 5 = Co = C2 = 0, i.e. 

i?:,R(/3, K V) < Er,R{/3, h, V) < El^^{f3, h, V) + C^h-\\\V\ff^l + r'^). (2.9) 

Note that norms of V in this Lemma are in the whole space. 
Proof Let H = Hh{A) =Th{A)-V. Since 

Tr [ilJrHilJr]_ > Tr [llJr[H]-tlJr]_ = Tr Vr[^]-^r, 

the first bound in (12.81) is trivial. 

In order to prove the second inequality, we write 7 := l(_oo,o](-f^) and calculate, 

Tr = iTr Oj^^H]. + [H].^',) = ^Tr ((^^^ + H^',)j) 

= ^Tr {{[H,iPr],^Pr]+2^lJrH^^Jrh). (2.10) 
Therefore, by the variational principle, and since [[if, ■i/'r.], "^r] = —2h'^(Vipr)'^ 

Tr [i/jrHipr]^ < Tr ipr[H]^ + h^Ti {Viprf-f. (2.11) 
In order to estimate the last term we apply a Lieb-Thirring inequality from |LLSj : 

Theorem 2.4. ILLSf There exist a universal constant C such that for the semiclassical Pauli 
operator Th{A) - V with a potential V G L^/'^{E?) fl /.''(M^) and magnetic field B = V x A e 
L^(M^) we have 

Tr [T,{A)-V]^>-Ch'' I [V]f-C[h''l \B?f\j [V]\)"' ■ (2.12) 

Using Tr 7/7 < and applying (I2.12p to the operator H — {ViprY, we obtain 
-Tr iVtl^rf^ >Tt {H- (V7/;,)2)7 (2.13) 



> 
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From 02.111) we thus have 



ErAP^ V) < inf [Tr t^^. [n{A) - V] ^ + {1 + 6)13 j \V x A 



+Ch-'(^ J [V]f + r-2) +C/i2r 3/3-3 J [V] 



4 1^-5 



Exactly as in the upper bound in Lemma [2m we can replace J^^ |V x with /gj-^-j |V ® 
at the expense of increasing /3 by a factor. This proves Lemma [2.31 for the Pauli case. For the 
Schrodinger case we can use the non-magnetic Lieb-Thirring inequality, i.e. the second line of 
f l2.13p simplifies to 

-Tr {VAYl > -Ch-' j [V+ (V^.)'] f 
and the rest of the argument is the same. □ 



3 Refined local semiclassics with scales 

We consider a local version of Theorem 11.41 for a model problem living in the ball B{i) of 
radius £ > 0. Note that not only the Hamiltonian but also the magnetic field has been 
localized, albeit on a twice bigger ball. We introduce two scaling parameters, i and /, that 
describe the lengthscale and the strength of the potential and we follow the dependence of the 
error term on these parameters. This scale invariant formulation will be convenient for the 
multiscale analysis in |EFS3] and also in the reduction of the following theorem to its special 
case. Theorem 14.11 Readers who are interested only in the semiclassical aspect of the result 
can think of / = £ = 1. 

Theorem 3.1 (Semiclassical asymptotics). There exist universal constants rzo € N and e > 

such that the following is satisfied. Let kq, /, i,ho > and let k, < KQf~^i~^ . Let ip G C^(M3) 
with suppip C B{i) and let V G C°°{B{i)) be a real valued potential satisfying 

< c„ri"i, < c„f ri"i (3.1) 

for every multiindex n with \n\ < uq. Then for the Pauli operator Hh{A) := Th{A) — V, 

inf (Tr [ijHh{A)ij]^ + 4? / I ^ ® ^1') - 2(27r/i)-3 f ^[qf [p^ _ V{q)\ dqdp 

< Ch-^+' f^-^^-' (3.2) 

for any h < hofi. The constant C depends only on kq, Hq and on the constants C^, in ( 13. ip . 
The factor 2 in front of the semiclassical term accounts for the spin and it is present only for 
the Pauli case. 
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Remark 1. For convenience, we introduced a new parameter k = I3~^h~'^ instead of /3. In 
the limit we consider, h — )■ oo, fSh"^ > const > 0, the new parameter k will remain bounded 
uniformly. We note that the constant factor /3 in front of the magnetic energy does not 
necessarily have to scale as h~'^ (see |EFS1] for other scalings) but our choice here is dictated 
by the application to the Scott correction. Similar results may be proven with a more general 
coefficient /3 ~ /i""^ with a certain range of 7 and the exponent of the error term will be 
e = e{'-f). We will not pursue the most general result in this paper. 

Remark 2. By variation of k, we obtain from Theorem 13.11 the following estimate 

I V ® A\^ < Ch-^+'f-'f-', (3.3) 

{21) 

for (near) minimizing vector potentials A. 

The following result can be viewed as a partial converse to (13. 3 p as it estimates the semi- 
classical error in terms of the magnetic field. Note that the assumption in (13. 4p below is much 
weaker than (13. 3p . 

Theorem 3.2. Let the assumptions he as in Theorem \3.1\ and assume that A satisfies the 
hound 

[ \V®A\^<Ch~^fi\ (3.4) 
Then, with e from Theorem \3.1\ we have 

Ch-^ff/^ [ \V (g)A\^y^\ch~^fi 

^ Jb(21) ^ 

> Tr [ipHhiA)^P]_ - 2(27r/i)-3 f ^{qf [p^ _ ^^(g)] _dgdp 

> Ch~^+'f^-'f-' - Ch~^fH [ I V ® (3.5) 

JB{2e) 

where the constants may depend on Hq and kq and on the constant in (13. 4p . 
3.1 Proof of the Main Theorem 11.11 

Using the local semiclassical asymptotics Theorem 13.11 it is an easy localization argument to 
prove the global result for compactly supported potentials. 

Recall the choice of a smooth cutoff function ip from Section |2l Let / = 1 and i > 1 
sufficiently large so that iIj^^x) := ipi^x/tj = 1 on the support of V . 
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The upper bound in Theorem II .11 is obtained by setting A = and using the first inequahty 
in Lemma 12.31 (with R = oo) 

E{l3, h, V) < 2Tr [-h'^A - V]^ < 2Tr ^-h'^A - V]^^)^ < 2Tr [i)i{-h^A - V)i)t]_, 

where the traces are on L2(]R^) and the factor 2 accounts for the spin. The upper bound now 
follows from the main result of [SSj which is formulated for = 1, but it clearly extends to 
any value of t. 

Theorem 3.3 ([SSl Theorem 12]). Let d > 3 and tp E Co"^^(R'^) be supported in a ball B of 
radius 1 and V G C^{B) be a real function. Let H = —K^A — V , h > {i, acting on L'^iW'-). 
Then, 

Tr LHM.4^Hi^]~ - (2^^)"' / - y{(l)]~ dpdg < C/i-'^+^Z^ (3.6) 

The constant C > here depends only on d, \\il)\\cd+i and \\V\\ci. 

For the lower bound in Theorem 11.11 we use the IMS localization with a partition of unity 
ipj + ipj = 1 and we borrow a small fraction of the kinetic energy to control the localization 
error. We have 

Th{A) -V>{1- 8)ij,{Th{A) - V5)ij, + 6[Th{A) - 8~^hH,] (3.7) 

with 5 e (0,|), h := iyilJt)'^ + iyi)iiY and Vs = {I - 5y^V . Here we dropped the term 
ipi{Th{A) — Vsjipi that is positive since = on the support of ipi. The second term in 03.71) 
is estimated by the magnetic Lieb-Thirring inequality similarly to fl2.13p : 

Tr [Th{A) - 6~^h^h]_ > -ChH~^/^ - Cp-^'h^d^^ - (3 j \W x Ap, 

where the constants depend on V and L Using ph"^ > Kq and assuming h < 1, we can combine 
the two error terms to obtain 

E{f3, h, V)>{1- S)El^{f3, h, Vs) - ChH~^ > ^,,2£(/3/2, h, Vs) - Ch^6~' - Ch^K 

In the last step we used the second inequality from fl2.8p with R = oo and then the monotonic- 
ity of Ei^R < in R. The energy Ee,2eil3 /2, h,Vs) is estimated in (with k = 2/3-^/1-^) 
and we get 

>/r3xR3 

where we used again that ipi = 1 on the support of Vs. It is easy to see that 
/ [p^-Vs{q)]^dpdq= I [p" -V{q)\ dpdq + 0{6) 

and finally, choosing 5 = h^^^ we obtain an error term Ch""^^^ assuming e < 1/4. This 
completes the proof of Theorem 11.11 □ 
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4 Refined local semiclassics with lower bound 



The main technical result of this paper is the following version of Theorem 13.11 where addi- 
tionally a lower bound on V, namely ( 14 ■2p is assumed. Furthermore, we impose the condition 
that V ■ A = 0. 

Theorem 4.1 (Semiclassical asymptotics) . There exist universal constants rio G N and e > 
such that the following is satisfied. Let Kq, /, i,hQ > and let n < Kq/"^^"^. Let tp G C(j"(M'^) 
with suppip C B{i) and let V G C°°{B{i)) be a real valued potential satisfying 

< c,ri"i, id^v] < c„/2ri"i (4.1) 

for every multiindex n with \n\ < no, and for some Cq > let 



mf{V{x) : X e B{i)} > Cq/^. (4.2) 



Then 



M'Jtt [t/jHf,{A)i>]_ + ^ [ |V® Ap) -2(27r/i)-3 f ^lj{qy[p^-V{q)]_dqdp 

<Ch-^+'f^-'f-' (4.3) 



for any h < hofi. Here inf^ denotes that the infimum is taken over all A G H^{M?) satisfying 
that 

V-A = 0, onB{M/A). (4.4) 

The constant C depends only on kq, Hq and on the constants C„, and cq in fl4.ip and fl4.2p . 
The factor 2 in front of the semiclassical term accounts for the spin and it is present only for 
the Pauli case. 

Note that this theorem is formulated in a scale invariant way, so it will be sufficient to 
prove it for / = £ = 1. 

We first show that the condition V ■ A = is irrelevant for the statement of Theorem 14. H 
i.e. one gets the same result by dropping the condition V ■ A = 0. 



Theorem 4.2. Suppose that Theorem 4-1 holds, where the infimum is taken over all A with 



V ■ A = on B{5i/A). Then Theorem\4.1\ also holds where the infimum is taken over all A. 



Proof. Clearly the unrestricted infimum is smaller than the one with the condition V ■ A = 
imposed. We will prove the opposite inequality, namely that there exists C > such that 

inf'fTr [^pHh{A)^P]_ + [ |V ® A\A < inf (Tr [^pHh{A)^p]_ + ^ [ |V ® A| 



(4.5) 
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where inf'^ means that we take the restricted infimum over A with V ■ A = on i?(5£/4). 
Since the constant k, is arbitrary this imphes the resuh. 

For an arbitrary A G H^, we can add a constant to A in order to get, by a Poincare 
inequahty, that 

[ A^<C [ |V® (4.6) 

Jb{2£) J B(2l) 

Clearly, this additive constant does not change the energy, i.e. on the right hand side of (14. 5 p 
we can assume that A satifies (14. 6p . 

Choose a localization function x ^ C*^(i?(3£/2)), x = 1 on i?(5£/4). Define Ay. to satisfy 

Vx A^ = E^:=Vx (xA), V-A^ = 0. (4.7) 

This system defines A^ up to an additive constant that we will choose in order to be able 
to apply the Poincare inequality below. Finally, define a vector field A' hy A' := x^x' fhen 
clearly V ■ A' = on S(5£/4). 
Then, 

Tr [i;Hh{A)ij]^ = Tr [ijHH{xA)i^]- = Tr [ijHh{A^)ij]. = Tr [ijHh{A')ij]., (4.8) 

where the middle equality follows from the gauge equivalence of A^ and x^i ci-nd the other 
two follow from the fact that x = 1 on supp ip. Also, using the Poincare inequality twice, once 
for A (14. 6 p and once for A^, we have 

/ \V®A'\^<f \v®A^\^ + cr^ [ \Ay\^<c[ |V®A^p 

JB{2e) J B{2e) J B{21) J B(2l) 

<C [ |V®A^|' = c/ Bl<c{r^ f A^+ [ \VxA\A 

Jr3 Jr^ ^ Jb(21) Jb{21) ' 

<C [ |V® v4|2. (4.9) 

JB{2e) 

□ 

4.1 Multiscaling 

To prove that Theorem 13.11 follows from Theorem 14.21 requires to perform a multiscale decom- 
position around the sets where V is too small and it violates (14. 2p . We present the setup for 
multiscaling for general potentials that will also be applicable for resolution of Coulomb like 
singularities. 

We will assume that the potential has a multiscale structure. Intuitively, this means that 
there exist two scaling functions, /, £ : — ?■ such that for any u G M^, within the ball 
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Bu{i{u)) centered at u with radius i{u), the size of V is of order /^(m) and V varies on scale 
i{u). Moreover, we also require that the continuous family of balls Bu{i{u)) supports a regular 
partition of unity. The following lemma states this condition precisely. This statement was 
proved in Theorem 22 of |SS] with an explicit construction0. 

We will use the notation Bx{r) for the ball of radius r and center at x and if a; = 0, we 
use B{r) = Bq^t). 

Lemma 4.3. Fix a cutoff function tp G C^(M^) supported in the unit ball B{1) satisfying 
J i/j"^ = 1. Let £ : — )■ (0, 1] be a function with ||V£||oo < 1- Let J{x,u) be the Jacobian 
of the map u [x — u)/i{u) and we define 

Then, for all x eM^^, 

[ M^f^iuy'^du = 1, (4.10) 
and for all multi-indices n G we have 

||9>„||oo < |n|=ni+n2 + n3, (4.11) 
where Cn depends on the derivatives of ip but is independent of u. 
We will require that the potential satisfies 

\d''V{u)\<CJ{ufl{u)-\''\ (4.12) 

for all n G uniformly in u in some domain C M^. In applications, Vt will exclude an 
/i-neighborhood of the core of the Coulomb potentials. 
For brevity, we will often use = l{u) and = f{u). 

4.2 Proof of Theorem [331 

We now show how Theorem 13.11 follows from Theorem 14. 1[ 

Proof of Theorem IJ.il By Theorem 14.21 we may forget about the condition ( 14.41) in Theo- 
rem |1]T] and consider the infimum over all A G H^. Moreover, since (13.21) is local, we may 
assume without loss of generality that V has compact support contained in i?(5£/4). After 
these remarks. Theorem 13.11 would follow from Theorem 14.11 once we have shown that the 
additional condition (14. 2 p in Theorem 14.11 can be removed. 

Since the statement of Theorem 13.11 is scale invariant we may also assume that / = £ = 1 
in the statement of that theorem. We will apply Theorem 14.11 and the partition of unity 
Lemma 14.31 with a different choice of / and i. 

^Multiscaliiig was introduced in semiclassical problems in [IS] (see also |Sob| ) 
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The statement of Theorem 14. II consists of a lower and an upper bound on the total quantum 
energy. Since the proof of the upper bound (see flS.ip ) does not use the condition fl4.2p . the 
upper bound in Theorem 13.11 follows immediately, so we only need to prove the lower bound. 
Define the smooth functions 

i(x) = f{xf:=^{V{xf + h"^Y^' (4.13) 

with some positive exponent with 2/5 < a < 1/2 and recall the notation = i{x) and 
fx = fi^)- Here K > (independent of h, but depending on V and ho given in Theorem 13. ip 
is chosen so large that 

||V£||oc < 1/4, 4<l/4, ueB{2). (4.14) 

By Lemma U3] there is a partition of unity ipu associated to i. Inserting this partition of unity, 
we have for sufficiently small values of h that 

Tr [i^{Th{A) - V)i^] _ + ^ [ I V ® A\^dx (4.15) 

Jb(2) 



= Tr 
> Tr 



+ ^ / IV® A| 



B{2) 



H?/ 

^P^[T,{A)-V-Ch'C]ij^ 

B(3/2) «n 



IV® A 



2 



B{2) 



> [ hi [MTh{A) -V- ChH-J)^^] +^ [ I V ® A\Mx} 

JueB{3/2) 1^ JBu{2iu) ' 



with '■= ipui^ and k' a fixed constant times k. We also used Tr [J 0„dn]_ > / Tr [0„]_dn 
for any continuous family of operators Ou- In arriving at f l4.15p we reallocated the localization 
error by using 

/ 3/ 

^<4<^, for all X G 5,(24) (4.16) 

which follows from ||V£||oo < 1/4. We also reallocated the magnetic energy by using the 
estimate 

|V® A(x)pdx5 < / / |V® v4(x)|2(24/3)-Mnda; 



ueB{3/2) JBui2eu) Jx&B(2) Ju€B^(4e^) 

= C [ \V®A{x)\Mx, (4.17) 

Jx(^B{2) 

where we again used f l4.16p . 

Consider now u e 5(3/2) and suppose first that |V^(ii)| > h"'- Then, if K in the definition 
of /, i fl4.13p is chosen sufficiently large, we have by the uniform bounds on the derivatives 
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of V that \V{x)\ > \V{u)\/2 > fl/A for all x G So ^ is satisfied (with cq = 1/4). 

Using (14.111) . the estimates in (14.1 p are also easily seen to be satisfied for the cutoff function 
il)u and for the potential V + Ch'^i~^ with £ = iu and f = fu (here we use that a < 2/3 and 
that K may depend on ho). We can therefore conclude by Theorem 14.21 that for any A, 

[ hr [UUA) -V- ChX')i^u] + 4^ / IV ® A\'dx]^ 

> [ / / MQYip' - v{q) - c/.^c^]_dgdj, - c/z-^+^e^^/n^ 

J{ueB{3/2y.\viu)\>hc} ^ {2'Khy J J ^ J ii 

> 2(27r/i)-3 / ^{qY [p^ - V{q)] _dqdp - Ch^-^'^h-'' - Ch-^+' [ C^'dw, (4.18) 

iR3xM3 Jb{3/2) 

where we used (14.101) and that h'^i"'^ < C/i^"^" to get the last inequality. Since i is bounded 
and a < 1/2, both error terms in (I4.18P are acceptable. 

For the set of u's where V{u) is small, i.e., |^(?/)| < h°', we use the magnetic Lieb-Thirring 
inequality. We introduce A' = x{A — c), where x ^ C^{Bu{2iu)), x = ^ on Bu{iu) and 
|Vx| < C/iu- Here c = (2e Then we have, using the Poincare inequality, 

/"|V®v4f<c/ (c'(^-c)' + |V® A]') <C" / |V®v4|2. 

So we get 

Tr [z^.(T,(A) -V- Ch'i-^)^^]_ + [ |V ® 

JBu(2iu) 

> Tr [T,(A') - (y + + 4^ / 1^ ® ^1' 



> -C 



[ [v + chX']f + ( [ [V + chX'YYUh-' ! \v®aA"'] 

^ JB^(iu) ^JBuitu) ' ^ JBui2£u) ' ^ 

1 



> -c{.'){h-' I [v+ chX']T + [ [y+ chx%] 

> -C(K')/i"'+'"/'^L (4.19) 

since \V{u)\ < /i" and 4 ~ /i", so we get \V{x)\ < \V{u)\ + < Ch"" for any x G Bu{D- 
In summary, we have 

/ {Tr [V^„(T,(A)-y-C/.X~')V^«]_ + 4^ / iV^Apj^ 

J{ueB{3/2):\V(u)\<h'^}^ l^'T' J Bu{2lu) ^^u 

> -C7/i"^+^"/2. (4.20) 
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Since a > 2/5, this error term is acceptable, in fact a = 4/9 optimizes the errors from fl4.18p 
and fl4.20l) . which shows that e < 1/9. Thus the final e in Theorem 13. II is the minimum of 1/9 
and the e obtained in Theorem I4.1[ This finishes the proof of the lower bound and therefore 
the proof of Theorem I3.1[ □ 



4.3 Proof of Theorem [372] 

We now show how Theorem 13.21 follows from Theorem I3.1[ 

Proof of Theorem \3.2[ By the argument given in the proof of Theorem 14.21 we may assume 
that V ■ A = on a neighborhood of supp^/^. Also, by a rescaling we may assume that 
f = i=l. 

The lower bound in (13. 5 p follows directly from Theorem 13.11 The upper bound in (13. 5p 
will be constructed using an explicit trial state. Define 



7 := 



f x«exp(-ic„x)7„exp(-2c„a;)x«7^ (4.21' 



with a length scale Lq < 1 that will be optimized at the end of the proof. Here Xu{x) = 
x({x — u)/Lq), with X ^ C^{B{1)) a real, positive function with J = 1- Clearly Xu is 
supported in the ball Bu{Lq). We choose the parameter c„ := ^^^^ A. Also 

In = l(-oo,0] {Xui'HhiA = 0)lpXu)- 

Notice that '-/u is real and acts like a scalar in spinor space. Therefore the contribution of the 
cross terms [cr ■ (— z/iV)][cr ■ A] + [cr ■ A] [cr ■ (— z/iV)] on ipXuluXuil' vanishes and we get 

Tr [xui^Hh{A - Cu)iJXu7u]Y^ 
= / (Tr [xu^Hn{A = 0)fe]- + Tr [xli^'{A - c^^j) 

Jb{2) ^ ^ -^0 

The du integration can be restricted to -8(2) by the support properties of V' and Xu and by 
Lq < 1. By a rescaling to a ball of unit size and an application of an optimal semiclassical 
result (recalled below as Theorem 15.60 we get that 

Tr [xu^H,{A = 0)i;Xu]~<2i2iTh)-' [ xl{qmQf[p'-Viq)]_dqdp + C{h/Lo)-\ 

(4.22) 

with a constant C independent of u. Performing the integration over u we will therefore get 
the desired upper bound if we can estimate the A^ term as well. For this, we apply Holder, 
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Lieb-Thirring and Sobolev inequalities to get (with 7^ := Xu'^lu'^Xu and p„ the associated 
density that is supported on Bu{Lq)), 



J Bu{Lq) J 



3/5 



where we defined D := —ihW for brevity. Here we used the choice of c„ and the Poincare 
inequahty in Bu{Lq) to control A — in by //oil V ® A\\i2(^Bu{Lo))- 
Since we have the following bound, which we will prove below, 

Tr [D^7j < Ch^^Ll (4.23) 

we can estimate 

/ Tr [xl^\A - c„)^J^ < CLlh-' / || V ® Al|i.(^^(,„„^ 

< CLlh-^ I |V® A|2. (4.24) 

Jb{2) 

Now we choose Lq by optimizing the error terms CL^^{h/Lo)-^ from (02]) and from (jlJl]). 
Taking into account Lq < 1, the result is the choice 

Lo = mmh, h^/^( f |V®A|2j (4.25) 

Therefore, we obtain the desired upper bound in (13. 5p . 

It only remains to establish fl4.23p . Applying (I4.22p for an upper bound (and using that 
the first term on the right hand side is non-positive) and the non-magnetic Lieb-Thirring 
estimate for a lower bound we get 

CLoh-' > Tr [xu^HhiA = 0)^xj- 

> -Tr [D^%] - Ch-'^ I y^/l ^4 26) 

Since L^^ < Ch~^ from (I4.25p and (13. 4p . this estimate clearly implies ( I4.23p . This completes 
the proof of Theorem 13. 2[ □ 
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5 Proof of the Semiclassical Theorem 



an 



The rest of the paper is devoted to the proof of the main technical result, Theorem 14.11 We 
can assume that f = i = 1 since f l4.3p is scale invariant. 
For the proof of the upper bound, we just set A = 0, i.e. 



where the last estimate is a consequence of Theorem 13.31 Actually, the exponent in the error 
term in the upper bound in (15. ip can be improved from —2 + 1/5 to —1 by appealing to 
Theorem 15.31 below. However, we will not need this. 

The proof of the lower bound in Theorem 14.11 has several steps. First, in Section [5TT| we 
localize onto balls of radius Lq ^ h^^"^ at the expense of an acceptable error of order 0{h^~^^), 
in fact we will choose Lq := /i^/^""^" with some small Eq, and e will depend on Eq. (Note that 
here Lq is chosen differently than in the proof of Theorem 13.21 ) Then we wish to replace A 
by a smoothed out version Ar on scale r := h^/'^^P <C /i^/^ for some small p > with an error 
of order 0{h^^'^), where e will depend on q. This will eventually be achieved in Section [5^ 
(Theorem 15. 2p . In order to do that, we will need apriori bounds on the magnetic field and 
the momentum distribution of the low energy trial states. To obtain these apriori bounds, in 
Section [6l2] we will introduce a second localization on an intermediate scale Li = /i^/^+eo ^j^h 
r ^ Li ^ Lq, and we show that on this scale A can be neglected. These apriori bounds will 
have a weaker precision of order 0{h^~^) since localization onto a short scale Li <^ h^^"^ is 
expensive, but this will be sufficient as an input for Section 15.21 

Finally we go back to the larger scale Lq and use a semiclassical result for the operator with 
the smoothed vector potential Ar. Note that the scale of regularity, r, of the vector potential 
is much smaller than Lq, so it is not a straightforward application of standard semiclassical 
results with C°° data. However, one can keep track of the dependence of the error terms in 
the standard semiclassical statements on the derivatives of the symbol, which in our case will 
be powers of Lq/t. If L^/r is not too large, this error can be compensated by the smallness 
of A, still rendering a few derivatives of Ar bounded which is sufficient for the semiclassical 
asymptotics. 

5.1 Localization onto balls of size Lq 

We introduce a partition of unity on the lengthscale Lq = < 1/4 with some sufficiently 

small £o > 0, i.e. we choose 0* G C^(M^), / (pl = 1, supp 0* C -B(l), and define 



inf'(Tr [^Aff,(A)^]_ + -^ 




3(2) 





'0 
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then 



-^0 



Inserting this identity into Tr [ipHh{A)ip]^ and using IMS localization, we have 

du 



Tr [iljHh{A)ij]_ = Tr 



B{3/2) -^0 



buHh{A)(l)u - h'\V(l)u\ 



> Tr 



du 



B(3/2) 



--^^<Pu{Hh{A)-Ch''Lf)<Pu^ 



> 



du^ 
— Tr 

B{3/2) Ll 



hu{Hn{A)-Ch^Lf)<Pu^ 



after reallocating the localization error. Notice that the du integration can be restricted to 
5(3/2), otherwise V'^ti = 0. We can again redefine the potential 

V :=¥ + Ch'^Lf. 

The change of the semiclassical term due to this modification, 

2{2Tih)-'^ j I ^{qf[p'^-V{q)\_dqdp-2{2Tih)-^j j ^{qf[p^-V^{q)\_dqdp 



< Ch"'h'LQ' < Ch~'+''\ (5.2) 



can be incorporated into the error term in (14. 3p . Therefore to prove the lower bound in (14. Sp . 
it is sufficient to prove that, for some e > 0, 

inf'( / ^Tr [,|J<P^Hi{A)<P^^|J]^ + 4^ / |V ® 

> 2(27r/i)-3 j j ^{qf [p^ - V+{q)\ dqdp - Ch-^+', (5.3) 

where := — V^. Since satisfies the same bound as V (second bound in (14. ip with 
/ = £ = 1), we can drop the upper index + and we will focus on the proof of (15. 3p . 

We reallocate the magnetic energy and consider the infimum over A for each u separately. 
Reallocation changes the coefficient of the field energy by a universal constant factor c > 
using the inequality 



B(2) 



\V®A\'>cf |V 

Jb(3/2) -^0 Jb^{2Lq) 



A'. 



Therefore, it is sufficient to prove that for each fixed u G i?(3/2) and c > we have 



Bu{2Lo) 



> 



2{2Txh)-' j j ^{qf(ti[p'-V{q)\^dqdp-Ch-^+^Ll (5.4) 
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and then integrating this inequahty w.r.t. cIm/Lq over u G 5(3/2), we will obtain f l5.3p . Recall 
that inf^ denotes infimum over vector potentials with V ■ A = on 5(5/4), in particular they 
are divergence free in a neighborhood of the support of tp. 

Defining (j) = i/j(j)u, Theorem 14. II will immediately follow from the following theorem 

Theorem 5.1. Let Lq = with a sufficiently smaller > 0. Let be supported in B(Lq), 

with < CnL^^^y Let the potential V satisfy 

\d'^V\ < Cn on B{Lo), and V{0) > Cq (5.5) 

with some positive constants Cn and Cq. Then for any n > we have 

inf'lTr [0i7,,(A)0]_ + 4^ /" \V ® A\A 
^ \ i^h^ Jb{2Lo) J 

> 2{2iTh)~'^ J J (Piqf - V{q)] dqdp - C/i^+""Alo (5.6) 

with a constant C depending on n, C„ and cq, where we recall A^g = /i~^Lq. Here inf^ is 
taken for all vector potentials satisfying V ■ A = on i?(5Lo/4). 

In the next subsections of Section [5] we will prove Theorem 15.11 We first smooth out the 
magnetic vector potential. A, on a scale r -C h^^"^. Theorem 15.21 in Section 15^2) states that the 
error of this replacement is negligible. This is the main technical result and it will be proven 
in Section ini Once A is replaced with a smoothed version A,., and we estimate the size of its 
derivatives, we can apply a semiclassical result (Theorem 15. Sp to evaluate the r.h.s. of (15. lip . 
Combining Theorem 15.21 and Theorem 15.31 will then yield Theorem 15.11 

5.2 Replacing A with a smooth version on Lg-box 

To prove Theorem 15. H we first smooth out the magnetic potential on a scale r = h^/"^^^ <^ h^/"^. 
Using the lower bound from Theorem l3.31 it is sufficient to consider magnetic vector potentials 
A such that 

8{A) ■= Tr [(PHh{A)<P]^ + 4^ / |V ® < f (0). 

'^'^ Jb{2Lo) 

We can assume that A has zero average, i.e. 

/ ^ = 0, (5.7) 

JB(2Lo) 
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since S{A) = S{A — c) for any constant shift c G in the vector potentiaL To see this, note 
that for any cutoff function ip, by the variational principle, 

Tr mn{A) - = inf {Tr ^(T^(A) - V)ij^ : < 7 < 1} 

= inf {Tr tp{Th{A - c) - \/)^/;e-^'^-^/ V"'^'' : < 7 < l} 

= inf {Tr ipiThiA - c) - V)ip-f : < 7 < l} 

= TT[ij{Tj,{A-c)-V)^]_. (5.8) 

In fact, Tr [ip{Th{A) — V)ip]^ is invariant under any gauge transformation, A ^ A — \/<f, 
ip -.R^ ^R, but /5(2io) I V ® is not. 

Now we state the main result of this section: 

Theorem 5.2. Let k, > be given. Let Lq = h^^'^~^" for some sufficiently small Sq > 0. Let 
(f) with supp0 C B{Lq) with \d'^(j)\ < C„Lq Assume that the potential V satisfies 

\d'^V\<Cn onB{Lo) and V{0)>co (5.9) 

with some positive constants Cn and cq. Let A be a vector potential such that S{A) < S{0) 
and satisfying (15. 7p and V ■ A = on i?(5Lo/4). Then for a := 1 — Seq we have 

h-^ I |V ® < C/i"Aio . (5.10) 

Jb{2Lo) 

Moreover, for p > sufficiently small there exists Eq > (in the definition of Lq ) such that 
with r := h^^'^'^P we have 

Tr {(t>[Th{A) - V]0)_ + 4:2 / |V ® > Tr (0[T^(A,) - V]<P) _ - CH^^'^Kl,- (5.11) 

l^l^ J B{2Lo) 

Here Ar = A * Xr is a smoothed out version of A on the length scale r and the constants in 
the estimates depend on the fixed unsealed cutoff function x, on n, Cn and cq in (15. 9p . 

We remark that this theorem involves only the potential in B{Lq). However, under the 
conditions (15. 9 p one can extend V to B{2Lq) with similar bounds on the derivative. In the 
sequel we will thus assume that V is defined in B{2Lq) with 

\d''V\<Cn on5(2Lo) and V{0) > Cq, (5.12) 

for some constants C„, and cq that is determined by, but may differ from, the constants in 
(15. 9 p and with a slight abuse of notations we will continue to use the same letters. 
Theorem 15.21 will be proved in Section [6] below. 
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5.3 A precise semiclassical result 

We state the following simplified version of [11 Theorem 4.5.2] (see also |l2l Theorem 5.2.2]) 
for reference. Recall that D = —ihW. 

Theorem 5.3. Suppose that (p is a localization function with 

supp0c5(l), (5.13) 
and H = {D + — V is a self-adjoint magnetic Schrddinger operator acting on L^(M'^) with 

Id^'Vl < c, Id^'Al < c, and |a"(/>| < c (5.14) 
for all multi-indices n G N'^ with \n\ < K . Then, for all h G (0, 1], 



TV ^2(1,3)02 [(^D + Af-V]^- {27ihr^ / / 02(g) [{p + A{q)f - V{q)\ _ dgdp 



< Ch~\ 

(5.15) 

Here K is a universal constant and C only depends on the constant c in fl5.14p . 
The same result holds for the Pauli operator Th{A) = [a ■ {D + A)]^ as well 

Tr 02 [Th{A) -V]^- 2(27r/i)-3 j j (P\q) [{p + A{q)f - V{q)\ _ dqAp\ < Ch-\ (5.16) 

Remark 5.4. By a simple shift of variables p ^ p + A{q) for each fixed q, we obtain that 

i2nh)-' j I <p'iq)[{p + Aiq)f-Viq)]_dqdp=i2nh)-'J J <P'{q)[p' - Viq)] _ dqdp, 

in other words, the presence of the magnetic vector potential plays no role in the semiclassical 
formula. 

Remark 5.5. Theorem \5.3\ in Ivrii's book is formulated for Schrddinger operators, as stated 
in fl5.15p above. For the purpose of semiclassical analysis, the Pauli operator, written as 
Th{A) = (D + A)2 + her ■ B — V , can be considered as a Schrddinger operator with 2x2 
matrix valued potential that is subprincipal. Therefore the analysis of Ivrii goes through for 
the Pauli case as well and this gives (15.161) . Alternatively, one can apply Ivrii's result for the 
Dirac operator Theorem 7.3.14] or Jl^ Theorem 5.2.23] and use that in the large mass m 
limit, the square of the Dirac operator converges to the Pauli operator. Although not stated 
explicitly, the error estimates in the cited theorems are uniform as m oo. In this way, 
taking the semiclassical limit followed by the large mass limit, the semiclassical estimate for 
the Pauli operator can be deduced from Ivrii 's result on the Dirac operator. 

We also need a slight modification of Theorem 15.31 when the localization function is 
inside the negative part: 
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Theorem 5.6. Suppose that (p is a localization function with 

supp0c5(l), (5.17) 
and H = Th{A) — V is a self-adjoint Pauli operator acting or L^(R'^, C^) with 

< c, |(9"A| < c, and < c (5.18) 

for all multi-indices n G N'^ with \n\ < K . Then, for all h > 0, 

Tr [0(T;,(A)- V)0]_-2(27r/i)-=^y" J (f)\q)[p^ -Viq)]_dqdp\ < Ch-\ (5.19) 

Here K is a universal constant and C only depends on the constant c in fl5.18p . Similar 
statement holds for the magnetic Schrddinger operator. 

Proof. Since this estimate is local, we may assume that supp V,suppA C B(2). For a lower 
bound we estimate 

Tr [<P{n{A) - V)4>]^ > Tr [<P{n{A) - = Tr [0(T,(A) - V)_<p]. (5.20) 

In order to prove the upper bound, we write 7 := l(_oo,o](-^) and calculate, 

Tr <P'[H]^ = ^Tr {<p'[H]^ + [HW) = ^Tr {{<p'H + H<p')^) 

= lTr (([iJ,0],0] + 20iJ0)7). (5.21) 
Therefore, by the variational principle, 

Tr [(/>iJ0]_ < Tr + /i^Tr {V(f)f-f. (5.22) 

In order to estimate the last term we apply a Lieb-Thirring inequality to the operator H — 
(V0)^ Using Tr H-f < 0, this yields 

-Tr (V0)% >Tt {H- (V(/>)^)7 

>-Ch^'l {[V]^ + {Vct>ff"-C[h~^ j \V^A\^f\j{[VU^{Vct>fYf' 

> -Ch'^ - Ch-'^'^. (5.23) 

For Schrodinger a similar Lieb-Thirring inequality holds but without the magnetic term. 

Combining the upper and lower bounds the estimate fl5.19p follows from fl5.15p in case 
oi h < 1. Finally, for h > 1 (15.190 trivially holds since the semiclassical integral is of order 
which is smaller than Ch~^. The quantum energy can be bounded by the Lieb-Thirring 
inequality after using fl5.20p : 

> Tr [(p{ThiA) -V)(f)]_>TT {Th{A) - V) _ > -Ch-^ - Ch-^''^ 

for the Pauli case and the term is absent for the Schrodinger case. Both terms are 

smaller than the error bar Ch~^ in fl5.19p . □ 
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Now we explain how Theorem 15.31 can be apphed to prove Theorem 15.11 from Theorem 15.21 

To estimate the right hand side of (15. lip from below, we first move the localization outside 
the negative part by the simple inequality 

Tr {<P[{n{Ar) - > Tr (j)[n{Ar) - V]^<j) = Tr (j)^[n{Ar) - V]^. (5.24) 

By unitary scaling x — Lqx, we scale the unit ball to the ball B{Lq) and we have 

Tr (j)^[Th{Ar) - V]. = Tr ^^[n/Lo{Ar) - V]., (5.25) 

where 

0(a;) = 0(Lox), supp0c5(l), (5.26) 

and 

V{x) = V{Lox), Ar{x) = Ar{Lox). (5.27) 

Notice that the semiclassical parameter has changed from h to knew '■= h/L^ = /i^+'^o which 
is much smaller than h^^"^. Theorem 15.31 will be used for the data with tilde with hnew and it 
provides a precise semiclassical asymptotics with a relative error term of order Ch"^^^ compared 
with the main term. In terms of the original h, this error is of order C/i^"*"^^". 

We now check that the derivative estimates (15.141) hold for V, Ar and </>. Notice that we 
may replace V and Aj. by localized versions in the left hand side of (I5.24p . Therefore, we only 
need to control the (finitely many) derivatives in (I5.14p in the ball 5(3/2). 

Clearly, the derivatives of and V are bounded since Lq < 1. For Ar we have the following 
estimate which will be proved as a consequence of Lemma 16.161 in Section 16.31 

Lemma 5.7. Assume that /g(2Lo) ^ ~ ^> then the following estimate holds, for \n\ > 0, 

_ / r„\ 2max{|n|-l,0) /" 

|9"Af<C„(y) V/ |V®A|2. (5.28) 

Lemma [5.71 combined with f l5.10p . yield 

~ / T/n \ 2 max(|n| — 1,0) 

<C(^) Llh--\ (5.29) 

B{3/2) V r / 

In order to use Theorem 15.31 we need uniform bounds on K derivatives of Ar- By Sobolev 
inequalities this corresponds to K' > K + 3/2 L^-derivatives. Inserting |?t,| = K' and the 
definitions of Lo,r and a in (15.291) . we get 

/ I^M^p < c/ii-2(^'-i)p-(2i^'+2)eo^ (5.30) 



5(3/2) ^ ' JB{2Lo) 
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for all n G with \n\ < K' . Since K' is universal, the right side of fl5.30p is clearly bounded 
for sufficiently small p and Eq, i.e. if 

1 > (2Ji' - 2)p - (2Ji' + 2)£o. (5.31) 

By Theorem 15.31 applied to the Pauli operator, we get for such values of p, Eq that 



Tr (t^'in/LMr) - V]- = 2(27r(/i/Lo))-^ j j <P\q) [{p + A^f - V{q)\ _ dqdp+0{{h/ L,)-') 

= 2{2nhr' j I ^\q)[p^ - V{q)]_dqdp + 0{h'+'''^AL,). (5.32) 



Recalling that the error in (15. lip was 0{h^~^^"ALf^) this completes the proof of Theorem 15.1 



□ 



6 Proof of Theorem 15.2 



6.1 Reduction to a constant potential 



We use the notations and assumptions of Theorem 15.21 the extension of V given in (I5.12p . and 
assume that S{A) < S{0). From (I5.12p we can assume > for h sufficiently small. Let x be 
a nonnegative, smooth, symmetric, cutoff function on with j = 1 and suppx C -B(l). 
We choose a new lengthscale Li < ^Lq such that h < Li < h^/'^ \ \ogh\~'^. For any u e we 
denote 

Xu{x) = XuMi^) = ^(~^)' J J3^l(^) = 1' (6-1) 

and we will mostly omit the lengthscale Li. In the applications, Li = h^/'^'^^° . 

We first prove a rough lower bound on the left hand side of (15.111) . This will eventually 
be used to get apriori bounds when we prove (15. lip . Along the way, we also prove (I5.10p . 



Theorem 6.1. Fix k > and assume that the potential satisfies (15. 9p . Let G C^{B{Lq)) 
with \d"'(j)\ < C„Lq and let A be a vector potential satisfying (15.71) and 'V ■ A = on 
S(5Lo/4). For any a < 1 and S > there is a constant Cs^a such that if h < hs, then 



Tr 



0[T,(A) - y]0 



K 



IV® A| 



B(2Lo) 



> / Tr 



Xu<l>[D'-V{u)]^(j)Xu 



du 
Jif 



- (Cs,ah'' + Ch'Lf + CLiyLo- 



Jb(2Lo) 

(6.2) 
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Moreover, we also have 



inf < Tr 



Jb(2Lo) 



< / Tr 



dn 



Xu(t)[D^ -V{u)\_(t)Xu j^ + C{h^L^' + Ll)AL,. (6.3) 



LI 



Remark 6.2. We have explicitly, for any u G 



Tr 



157r' 



:h--'^[V{u)fl' / 



and thus 
After expansion 



Xu(t^[D^-V{ 



u 



dn 



157r2 



/ / [V{u)f'\l{x)<P{x)Mu. 



[V{u)f^ = [V{x)f^ + V[V{x)f\x -u) + 0{[x - 
and since, for any x G M^, 

j (^uxlix){x -u) = 

by symmetry, so we have 



Tr 



Xu<p[D'-V{u)]_<l)Xu 



du 



^-3 / ^5/202^^(^2^-3) / 0: 



Lf 157r2' 

Moreover, we have the following upper bound on the magnetic energy: 
Corollary 6.3. Under the conditions of Theorem \5.2\ we have 



2k} 



for all sufficiently small h. 



^1'' Jb{2Lo) ^ ^ 



(6.4) 



Corollary l6.3l is an immediate consequence of f l6.2p from Theorem 16.11 choosing 8 sufficiently 
small, using Ei^A) < S{0) and Theorem 15.61 (with rescaling it to the ball of size Lq). Choosing 
Li = /i^/2+'^o g^j^^ Q, = 1 — 35^^ we immediately obtain f IS.lOp . the first claim of Theorem 15.21 



Proof of Theorem \6.1\ First we give the lower bound (16. 2p . Recall that we can extend the 
potential V from B{Lq) to B{2Lq) by keeping similar derivative bounds as in (15. 9p . see (I5.12p . 
By the IMS formula, and Taylor expansion 



V{x) - V{u) - VV{x) -{x-u) < C{x - u) 



(6.5) 
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for all X G supp 0, m G supp (p + Li supp x, we have 



0[r,(A) - r]0 = / x^(x)0(x)[r,(A) - v{xmx)^ 



du 



> I Xuixmx) [n{A) - V{x) - Ch'L^'^<f){x)xu{x)jj^ 



du 



> I Xuix)(l){x) [n{A) - V{u) - VV{x) -{x-u)- Ch'Li' - CLf J 



Xu{x)<p{x) [niA) - V{u) - Ch'Ll' - CLi\<P{x)xu{x)jj^ 
In the last step we used that 



du 
1 ' 



(z — u)x^{x)du = 0. 
Using that Tr [^^. > ^ . Tr we get 



Tr 



(j)[niA) - v](i> 



> / Tr 



Xu^i^) - V{u) - Ch'L-,' - CL'Mxu 



du 



(6.6) 



(6.7) 



(6.8) 



In Theorem 16.41 of the next section, we will prove that A can be neglected. In order to 
facilitate the estimate, it is convenient to ensure that A has zero average on the ball By_(2Li). 
We define 

A. (6.9) 

Bu{2Li) 

Similarly to the argument in (15. Sp . we can subtract c„ from A in (16. 8p and we have 



Tr 



<p[n{A) - v]<p 



> / Tr 



Xu(l)[ThiA-Cu)-V{u)-Ch^Lf-CLj](t>Xu 73. (6.10) 



The M-integration in (I6.10p can be restricted to m G i?(|Lo) by the support properties of 
and Xu- For each fixed u G B{^Lq) we will use Theorem 16.41 proven in the next section to 



estimate the trace in the integrand. Define 

:= V{u) + Ch^Lf + CLj. 



(6.11) 



By fl5.12p we know that Co/2 < K < C. With the choice of L = Li, W = Vu and 77 = Xu4', 
we see from Theorem 16.41 that for any a < 1, e > and for any 6 > there is a constant Cs^e 
such that 



Tr 



Xu<P[Th{A - cj - V{u) - Ch'L^^ - CLl]<f)Xu _ 

> Tr Xu^[D' - Vu]_(PXu - SV^'^h-^ [ |V ® A\^ - Q./i^K'^^Al, (6.12) 

JBu(2Li) 
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as long as 
and 



|V ® < CKl.V^ 



(6.13) 
(6.14) 



Bu(2Li) 



hold for some constant C. The constant Cse will also depend on the constants in f l6.13p and 

(EH- 

Recalling Li = h^/'^^^°, with the choice a = 1 — 35o and e < £0/2, we see that fl6.13p 
is always satisfied since co/2 < Vu < C uniformly in m G B{2Lq). To guarantee the second 
condition f l6.14p for the availability of the estimate (I6.12p . we split the integral on the r.h.s. 
of (16.1 op as follows 



(6.15) 



'B{5Lo/4) 

where we defined 



6h"' 


1 


|v 




















|v 








JBu{2Li) 









The estimate (I6.12p will thus be available for u's in the first integral and it yields 

r "I (jlu 

l{u e S<}Tr [xu(l>[ThiA - c„) - Viu) - Ch^Lf - CL?]0x«J 

> / l{ueE^}(TTXu<P[D'-V^]_(t)Xu)^-C5h''' [ |V® A|2-C5>°Ai,. 

J ^ ^ -^1 Jb{2Lo) 



(6.16) 



With an explicit calculation, for any /3 > constant and cutoff function rj, we have 

2 



Tr 7][D' - I3]^r] = 2(27r)~'' / 7]\x)dx / [{hpY - /3]_dp 



ISvr' 



V 



(the additional factor 2 comes from the spin degeneracy). Since |V^(?i)| < C, /i < Li < 1 and 
> 0, we obtain 

K'/' < [V{u)]T + C{h'Lf + L?), (6.17) 

and thus 



Tr Xu(l>[D^ - K]_0X« > Tr Xu(I>[D' - V(w)]„0X. - Cih-'Lf + L2)A^^. 
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Inserting this into f l6.16p . we have the following estimate for the first integral on the r.h.s. of 

I l{nGS<}[...]^> I l{ueE^}(TrXu</>[D'-Viu)]_<Pxu)^^ 

- C5h~^ [ |V ® Af - [Cs,ah" + Ch^Lf + CLI]Al,. (6.18) 

Jb{2Lo) 

The traces in the second integral of fl6.15p are directly estimated using a Lieb-Thirring 
bound. For this estimate we will have to localize the magnetic field. Fix u G H> and choose a 
cutoff function (j)' with supp0' C i?„(3Li/2), |9"</)'| < CnL^"' and such that 0' = 1 on B{Li) 
and define, for the purpose of this proof, 



A' = {A~ c„)0', V'{x) = Vu ■ l{x e 



(6.19) 



Clearly 



/ |V®Af < / \W®A\'' + CLf! \A-Cu\''<c[ |V®A|2 (6.20) 

where in the second step we used the Poincare inequality. 

By the magnetic Lieb-Thirring inequality. Theorem 12. 4^ and f l6.20p : 



Tr 



Xu<P[ThiA - cu) - Viu) - Ch^L^^ - CLl]<Pxu\ _ = Tr 

> -ch~^ [ K'/' -c(h-' / |v X AfY'f [ v: 



x.0[T,(A') - v'^x^. 

1/4 



(Li) 
B„{2Li) 



(6.21) 



Bu{2Li) 



using that Vu < C and h < hg. In the last step we used m G S> to estimate V^Al^ by the 
local magnetic field energy. Integrating out this inequality for u G B{5Lq/4), we obtain the 
following bound for the second integral on the r.h.s in f l6.15p 



/ l{ueE^}[...]^>-C6h-' [ |V®A 

J ^1 JB{2Lo) 



(6.22) 



Finally, the missing piece of the non-magnetic term on the r.h.s of fl6.18p for u G H> can 
be estimated by the standard Lieb-Thirring similarly to fl6.2ip and f l6.22p 



1(m G S>)Tr Xu(p[D^ - K]0X« 
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IV ® Ar. 



(6.23) 



Su(2Lo) 



The estimates fl6.18p . fl6.22p and fl6.23p inserted into fl6.15p and fl6.1Up complete the proof of 

For the proof of (16.31) . we define the spectral projection 

In := l(-oo,o] {D^ - V{u) + C{h^Lf + Li)) 

and 



1 ■ j Xu'JuXu j^^- 

Note that < 7 < 1 by (16. ip . We will also use that the density of 7^ is given by 



An 



h"^[V{u)-C{h^Lf + Ll 



2M3/2 



Then, by Taylor expanding V up to second order, similarly as in (16. 5p but using W{u) 
instead of W{x), we have 



inf < Tr 

A 



<Tr (t)[D^ - 1/]07 



IV® > < Tr 



Tr {Xu(t>[D'-V](l)Xulu) 



\ du 



J LI 



< y Tr (xu<P[D' - V{u) - VV{u) ■ {x - u) + CLjjcPxulu) ^ 

= y"Tr (^^^^[D'-Viu) + CLl]^Xu7u)^^- J J 0,^ix)VViu) ■ (x - u)(i)\x)xlix)^^dx 

< / Tr {xl<P'[D' - V{u) + Ch'Lf + CL?]7„) ^ 



Cih~ 




[V{u) - C{h'L-' + Ll)]fvV{u) ■ (x - u)q^\x)xl{x)^^dx 



< y Tr {xu<l>[D' - Viu) + Ch'Lf + ^ + CLjh-' j , 

< y"Tr (xu<P[D' - V{u)Uxu)^^+ C{h'Lf + 



(6.24) 



In the last but one step we used (16.70 and that 

[Viu) - Cih'Lf + LDY^^VViu) = [Vix) - Cih'Lf + Ll)Y_l\vix) + 0(Li) 

if |a; — n| < CLi. In the last step we used an estimate similar to (I6.17p . This proves ( 16. 3p 
and completes the proof of Theorem 16.11 assuming Theorem 16.41 to be proven in the next 
section. □ 
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6.2 Removing A 

In this section we estimate the effect of removing the vector potential from the operator 
H{A) = Th{A) — W with a constant potential > on a scale L. In the applications, L will 
be Li = /i^/2+£o Let Hq = H{A = 0) = D'^ -W and define the free projection 

P = l(„oo,o](i^o). (6.25) 

We also set Al := h'^L^. 

Theorem 6.4. Given < a < 1 and < e < |(1 — a), let L be a lengthscale that satisfies 

<L< Ch^/W \ogh\y^W"^/\ (6.26) 

Let rj be a smooth cutoff function supported on B{L) with |9"?7| < CnL~^ . Let A G if]^j,(M^) 
be a vector potential. We will assume the gauge to be such that V ■ A = on B{5L/4) and 

/ A = 0. (6.27) 

We also assume that 

f^~2y^i/2 f \^ ^A\^ <CW^/^Al. (6.28) 

Jb{2L) 

Then for any 6 > there exists Cs,e, depending also on the constants in (16.26P and (16.28^ . 
such that for any density matrix < 7 < 1 we have 



Tr H{A)ri-fri = Tr [Th{A) -W]r]-fr] > Tr r][Ho]-V-Cs,eh''W^/^AL-C6h-^W^/^ [ \V®A\'^. 

J B{2L) 

(6.29) 



Remark 6.5. The gauge choice fl6.27p implies that one can use the Poincare inequality on 
B[2L) to conclude that 

A^ < CL^ [ |V® Ap. (6.30) 

B{2L) J B(2L) 



Proof of Theorem \ 6.4\ We start with localizing the vector potential which will be used 
later. Let (f)' be a standard localization function on B{5L/4), i.e. = 1 on B{5L/4), 
supp0' C B{3L/2), and define 

A' := (j)'A, (6.31) 

in particular V ■ A' = on B{5L/A). Note that this definition of A' is different from fl6.19p 
used in the proof of Theorem 16.11 in Section 16. Ij the prime notation will always indicate a 
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trivial cutoff outside of the appropriate domain we actually work on. In Section 16.21 we will 
use (Km . 

We have rjA = r]A', so the sole purpose of this modification is to guarantee that only the 
local part of A will be taken into account. The Poincare inequality f l6.30p remains valid in the 
form 

'[AfKl A^<CL^I \V(S)A\''. (6.32) 



JB{2L) JB{2L) 

6.2.1 Decomposition in energy space 

We introduce a dyadic decomposition around the non-magnetic Fermi surface using cutoff 
functions in energy space. The final result of this section is given in Lemma [6.71 

Let i > 0, he smooth cutoff functions such that suppxi C [| • 2*~^, | ■ 2*], Xi ^ I5 



X^it) = 1 for t G [f ■ 2-1, 1 ■ 2% \Vx^\ < C ■ 2- and 



Define cutoff functions on M'^ by 

/.(w) = X.(— fT^j for ^>0, ueR' 
\ wW / 

and 

fiiu) = X\z\ 777- for ^ < 0' ^ e ^^ 

V wW / 

with some w such that h < w < 1. Setting Iq := [| loggU^I] + 1, where [■] denotes the integer 
part, we clearly have /j = if z < — zq- Define Wi := 2^^^w, then /j is supported in a spherical 
shell of volume CwiW^^'^ for \i\ < iq. 
Clearly 

i>0 i<0 

SO we can define 

/O M = XO 777- 

V wW 

with an appropriate cutoff function Xo^ with < Xo ^ 1; |Vxo| < C , so that 



i.e., /o(m) is supported in the regime where \v? — W\ < ^wW and /o(m) = 1 where — W\ < 

fjj=0, if \t-j\>2. 



IwW. Note that 
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We also define /> by 

i>io 

We note that the support of />(t) lies entirely in the regime |t| > 2 and 

|V/>|<C, |supp(V/>)|<a (6.33) 
For each i with < |i| < io we also define enlarged cutoff functions fi by 



/>) = X.(^^) for ^>0 



and ^ 

fi{u)=X\i\{ TTT-) * < 0' 

V wW ) 



where Xi, i > 0, are cutoff functions such that suppxi C [| ■ 2* ^, | ■ 2*], Xj < 1, Xi(^) = 1 for 

' 4 



t e [f ■ 2^-1, f ■ 2'] and \Vxi\<C- 2-\ We also set 



V wW 

where Xo < 1) Xo(^) = 1 for \t\ < 2, suppxo C [—3,3] and |Vxo| < C. We can similarly define 
7>by 

i>io 

Note that x* = 1 on the support of Xii therefore we have 

hl = h 0<K|<^o, and fj> = f>. (6.34) 
These extended functions clearly satisfy 

fi<fi-i + fi + fi+i for \i\<io, (6.35) 
fio < /.o-i + h + /> (6.36) 

and < < 1, < 7> < 1. 

Finally we define the momentum cutoff operators by 

Fr.= MD), F^:=MD), Fr.= UD), F> := 7>(D), (6.37) 

then all inequalities fl6.34p . f l6.35p and f l6.36p are clearly satisfied as operator inequalities if 
the functions / are replaced with the operators F. 
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In the first step, we neglect the positive A"^ term: 

Tr H[A)'r]'-fri > Tr H[A)r)^r) 

with 

H{A) := Ho + (r{D)a{A) + (r{A)a{D) = Hq + 2D A + ha{B) (6.38) 

using V ■ A = on the support of rj and V x A = i? in the last identity and introducing 
the shorter notation (j{v) = cr ■ v ioi any vector v. Note that the formula f l6.38p holds if the 
kinetic energy Th{A) is the Pauli operator, for the Schrodinger case we would have the simpler 
formula 

H{A) := Hq + DA + AD = Ho + 2D A. 

Here we adopted the convention that DA = Xl^=i(~'^^^i)^j5 — Yl^j=i^j{~'^^^j) ^^'^ 
D ■ A = Yyj=i -ih{VjAj) = -ihdivA, in particular, [D, A] = D ■ A. 

Remark 6.6. We remark that neglecting A? is affordable since 7 is close to the projection P, 
defined in fl6.25p . and the density of the projection Qp{x) is hounded by Ch~^W^^'^ , thus 

Tr A^r]jr] ^ ! A^Qp = ch~^W^'^ [ A^ < cH-^LV^ [ | V ® 

J J Jb{2L) 

= {h-'L^W)h-^W'/^ [ I V ® A\' « h~^W'/^ [ I V ® A\^ 

J B{2L) Jb(2L) 

by the Poincare inequality fl6.30p and L <^ h^^'^W"^^'^ (see (16.261) ). 

By the reality of the projection P, and since it acts like a scalar in spinor space, we have 
Tr [a{D)a{A) + a{A)a{D) + a(5)]r/Pr/ = Tr {DA + AD)r]Pr] = 0. (6.39) 

Thus we get 

Tr H{A)T]'jr] >Tr H{A)r]-fr] - Tr [a{D)(r{A) + a{A)a{D)]r]Pr] 

=Tr F^Hor]-fr] + Tr [a{D)a{A) + a{A)a{D)]7]{'y - P)7] + Tr (1 - F^)Ho7]Pr] 
+ 5^ Tr F^HiA)rji^ - 1)^ " E F^H{A)t^{P - l)r/ 

i<0 i<0 

+ ^ Tr F^H{A)r]-fr] -^^r F^H{A)r]Pr]. (6.40) 

j>0 i>0 

The combination of the first and third term on the r.h.s of (I6.40p gives the main term in 

dssni): 

Tr FoHoFoVlV + Tr Ho{l - F^)r]P7] >Tr [t^FoHoFov]^ + Tr [r/(l - F^Y/^Ho{l - F^y/^ri]^ 

>Tr v{[FoHoFo]^ + [(1 - F^Y/^Ho{l - Fo^)^/^]_)r/ 
=Tr r][Ho\-V- (6.41) 
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Notice that the sum on the RHS of fl6.4UI) is real though each individual term might not be. 
To symmetrize some of the terms we take the real part and we have proved 

Lemma 6.7. With the notations introduced above, we have 

Tr H{A)r]-fr] >Tr r][Ho]_r] + 3ftTr F^[a{D)a{A) + a{A)a{D)]r]{'y - P)r] 

+ J2 F^H{A)7]{'y - - ^ 3?Tr F^H{A)r]{P - l)r/ 

i<0 i<0 

+ JJTr F^H{A)7]-f7] - J2 3f^Tr F^H{A)r]P7]. □ (6.42) 

i>0 i>0 

The main term is the first on the r.h.s of (16.421) . In the next subsections we estimate the 
other terms and we show that they can be included in the negative error terms in (16. 29^ . 

6.2.2 Estimate of the a{D)a{A) + a{A)a{D) term in (K^ 
Lemma 6.8. For any 6 > we have 

Tr F^ [aiD)a{A) + ^(A)a(D)] vil - P)v 

<C6-\h + wLW^'^ + h-WL^W]W'''^kL + ^h-^W^/^ f \V®A\^. (6.43) 

J B{2L) 

Proof. Using V ■ A = on B{5L/4) and t]A = rjA' (recall the definition of A' from (16.311) ) 
and the locality of the operator we have 

F^ [(y{D)a{A) + a{A)a{D)\ 7/(7 - P)r/ = 2F^DA'r]{-f - P)r] + hF^a{B')r]{^ - P)r] (6.44) 

with fi' = V X A'. Note that 



/ [Bf < [ A^+ [ 


[(j)f\V X A\^ <C 1 \V®A\^ 


(6.45) 


' J B{2L) Jr 


3 Jb{2L) 





using (16.271) and the Poincare inequality. 

We will estimate |Tr F^DA' r]'-f'r]\ and j/iTr FQa{B')ri^ri\ for any density matrix < 7 < 1, 
and then the same estimate can be applied to 7 replaced with P as well. We start with the 
second (magnetic) term in (16.441) . and in fact we prove the following stronger estimate that 
will be used later 

h\Ti FMB')v1V\ < E ^(Tr F^iBfy'ijT i^V^)'^' 

\i\<io \i\<io 

\i\<io 

< 5h-^W^'^ I |V ® A|2 + C6-^hW^'^XL. (6.46) 

Jb{2L) 
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Here we used 06.451) and that the diagonal element of Fi, \i\ < io, is bounded by 

supFi{x,x)<[ f,{hp)dp < Ch-'^WiW^'"^ (6.47) 

since the support of /j is a spherical shell of width CwiW^/"^ and radius of order W^^"^, where 
we recall Wi = 2'*'w. In particular, fl6.46p shows that the magnetic term in ( I6.44p can be 
estimated as required in Lemma [6.81 

For the first term in f l6.44p . with a general density matrix 7, we insert Fq and Go := 1 — -fo 
to get _ _ _ _ 

|Tr F^DA'r]-fr]\ < |Tr F^ D A' Forj-fT] Fq] + |Tr F^A'Gori-frjFoD], (6.48) 

where we used (16.340 and the cyclicity of the trace. We will estimate these two terms sepa- 
rately. 

The first term in fl6.48p is estimated as 

\Tt F^DA'FoVlvFol < (Tr [Af For/7r/Fo) (Xr Fo^Z^^^o^T^^o) 

/ _ \ 1/2 / N 1/2 

< C (Tr [AfF^j (wTt F^rf^ 

< 6^'h-WL^W'/^AL + 6h~^W'/^ [ I V ® A\\ (6.49) 

Jb{2L) 

Here we first used that is bounded by CW when multiplied by Fq. Then, similarly to 
fl6.47p , we estimated the diagonal element of Fq and Fq as 

supFo(x,x) < supFo(a;,a;) < Ch'^wW^^'^, (6.50) 

X X 

and finally we used fl6.32p . 

The second term in (16.481) is estimated as 

|Tr F^A'GqV1vFqD\ < (Tt F^ A' GIA' F^^ (Tt D^F^ir^^r^y^ . (6.51) 

For the second factor we can use the previous bound 

Tr D^F^{r]-fr]y < GWTi F^^f < Gh-'^wW^'^L^ = GwW^^^Al. (6.52) 
To estimate the first factor, we will need the following lemma, whose proof is postponed: 

Lemma 6.9. Let f and g be positive real functions on M.^ such that fg = 0. Then for any 
real function a G H^{M.^) we have 

Tr LHR4fiDH^)9iD)a{x)f{D)] < /i-'||/||oo||^7||oo|| V/||i||a||2|| V ® . (6.53) 
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Using that Fq = fo{D) with a function /q satisfying |V/o| < Cw |supp/o| < 

CwW^/^, thus llV/olli < CW, we obtain 

TiF^A'GlA'F^<Ch~^W\\A'\\2\\V®A'\\2<Ch~^Lw[ \V A\\ (6.54) 

Jb{2L) 

Here we used fl6.32p and the second inequahty in (16 .45 p . 

Combining these estimates and separating the two factors in (I6.5ip by a Schwarz inequahty 
we obtain the following bound 

|Tr F^A'Gori-friFoD] < C5~^wLW^Al + 5h"^W^^^ [ |V ® A\\ (6.55) 

Jb{2L) 

This completes the proof of Lemma 16.81 □ 
Proof of Lemma \6.9i By passing to Fourier space we get 



Tr L2^^s)[fiD)aix)giD)aix)f{D)]= // f{hp)a{p-q)gihq)aiq-p)fihp)dpdq. (6.56) 
We now use that fg = Oto rewrite the above integral as 

[/(M - fihq)]\a{p - q)\^g{hq)f{hp) dpAq 
< ||/||oo||^7||oo // l\h{q-p)-Vf{hp + th{q-p))\dt\a{p-q)\''dpdq 

\Vf{v)\dv [ IzWaiz)]"^ dz. (6.57) 



oo y oo 



The result now follows upon passing back in x-space in the integral over a. □ 
6.2.3 Error terms in (16.421) for i > io 

We now deal with the two sums in the last line of (I6.42p in the regime where i > io- 

Lemma 6.10. Under the conditions LW^^"^ <ti h^^^ and assuming (16. 28 p . we have for any 
fixed 6 > 



3ft VTr F^H{A)r]-fr] > -Ti [D^Flr]-fr]]-C6h-^W^/^ [ \V^A\'^-6Tt [D^Flr]-fr]] (6.58) 

^>io 2 >^^(2i') 

if h < hs- Moreover, for any N > 1 and h < kg, we also have 

\YTTF^H{A)r]P7] <C6h-^W^/^ [ l"^ ® + CnW''/^ ( . (6.59) 

Jb(2L) KLW^^^J 



1>'10 



' B{2L) 
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Proof. Using the first formula in fl6.38p . we write 
3? ^ Tr FfH{A)r]^r] = JJTr [Hq + a{D)a{A) + a{A)a{D)] r]-fr] (6.60) 

i>io 

> ^Tr [D^Flri-fT]] - |Tr F^a{D)a{A')r]-fr]\ - |Tr F^a{A')a{D)ri-fri\. 

Inserting G'> := 1 — F> into tlie second term, and using (16 .34^ . we liave 

|Tr Fla{D)a{A')7]-f7]\ < |Tr Fla{D)a{A')Fyr]-fr]Fy\ + |Tr Fla{A')G>7]-f7]F>a{D)\. (6.61) 
For tlie first term in (I6.6ip . with the notation a3> := FyT]'-fT]F^, we use 

|Tr Fla{D)a{A')Fyr]-fr]F:,\ < (Tt [Af w>y^^(^Tr D^w^^^l (6.62) 

Applying Holder, Sobolev and Lieb-Thirring inequalities and the bounds (I6.32p . (I6.45p . we 
have 

Tr[Af^^<{l[A']f\lfjf' 

KCL^/'f [ \V ® A\^) (h~^TT D^uS^\ (6.63) 

^ Jb{2L) ^ ^ ^ 

where ^>(x) := cD>(x, x) denotes the density of a;>. Thus (I6.62p can be estimated as 

4/5 



3/5 



\Tr Fla{D)a{A')Fy7]-f7]Fy\ < C L^'^^h-'^'^ ( j | V ® Al^V^Vxr 

^ Jb{2L) ' 

<Csih^^^^Lfh^^W'^ [ \V ® A\^ + 5Tt D^u> , (6.64) 

Jb{2L) 

where we used fl6.28p in the last step. Considering that LW^^"^ <^ h^^'^, we obtain that the 
first term in ( ]6.6ip is bounded by the two negative error terms in ( 16.58^ . 

For the second term in f l6.6ip . after a Schwarz inequality, we use Lemma [6.91 similarly to 

(EMD: 

\Tt Fla{A')Gy7]^7]Fya{D)\ < (Tt F^A'ClA'P^y^^ (Tt D^FHrj^r])'^^^^ 



< Cbhr^l'^W^I'^ [ I V ® + (5Tr D^F^rj-fr] (6.65) 

Jb(2L) 
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for h < hs, using LW^/'^ <^ h^^'^ and estimating ||V/>||i < CW . Therefore the second term 
in f l6.6ip is also bounded by the two negative error terms in f l6.58p . 
We now estimate the third term in (16.601) : 

|Tr Fla{A')a{D)r]-ir]\ < |Tr Fla{A')Fya{D)r]-ir]\ + |Tr Fla{A')G>a{D)r]-ir]\ (6.66) 

\V2/ \l/2 



< (Xr Fla^AfFlu^^ {Ti D^u^^ 
+ (Xr Fla{A')G>a{A')Fiy^ {Ti Fy7]^7]a{D)GMD)vivF>y^^ 

< c(Tr [AfFl^yFiy^^(TT D 



1/2 /_ xl/2 



C (Ti F^A'CyAFl^ (Tt 



1/2 /_ \l/2 



Here we used that F^ = F^Ft, and that a{D)G>a{D) = D^G> < D\ The second term in 
the r.h.s of (16.661) is estimated exactly as (16.651) . For the first term we essentially repeat the 
estimate (16.631) : 



Tr [AfF^UyFl < GL'^^ if | V ® aA (h-^Ti D 

^ Jb(2L) ' ^ 



3/5 



after estimating Tr D^F^uj^F^ < Tr D^Uy. This completes the proof of (16.581) . 

For the proof of (I6.59p . we write 

|Tr FlH{A)T]Pr]\ < |Tr F^HovPvl + 2|Tr F^ [a{D)a{A') + a{A')a{D)]r]Pr]\. 

Notice that in the estimate of the a{D)a{A') and a{A')a{D) terms on the r.h.s of ( 16.60p 
explained above is valid for any density matrix 7, in particular also for 7 = P, thus 



|Tr FlH{A)r]Pr]\ < (1 + 6)\Tt F^Hor]Pr]\ + G6h^^W^^^ [ |V ® A 

Jb{2L) 



(6.67) 



On the support of F> we have \Ho\ < 2D^, thus Tr F^Hor]Pr] < 2Tt DF>r]Pr]FyD. To 
estimate Tr DF^'qP'qF^D, we use (I6.72p from Lemma 16.111 below. We set i := hW^^^'^, 
f{p) ■= fy{W^/^p)p and g{p) := l(|p| < 1), so that Fy{D)D = W^/^f{-iiV) and P = 
gi-iiV) 

Tr DF^r)Pr)F^D = pF>r/P||L < CnW/'I^j^) (6.68) 

using that the supports of / and g are separated by a distance d of order one. The condition 
i < Ld is guaranteed by (I6.26p . Upon combining this with (I6.67P we have completed the proof 
of Lemma I6.1UI □ 



41 



Lemma 6.11. Let f and g be real functions so that their supports are separated by d, i.e. 

dist{supp{f), supp{g)) > d. (6.69) 

Let ri{x) := t]q{x/L) with L > 0, where rjo is a smooth function with compact support. Then 
for any N > and i < Ld we have the following bounds 



\\fHeV)vgHiV)\\Hs<C^{jj^) {j 



and 



l\f^_aV)vgHm\HS<CM{^) [-J 

fip) 



\\fHiV)vgHiV)\\Hs<C 



N 



dl 



1+p^ 



2||5'l|2 
siblloo, 

\g{p){l + py\ 



(6.70) 

(6.71) 
(6.72) 



where depends only on tjq and N. 
Proof. Since rjo G C^, we have 

|^o(p)|<C^(l+p')"^/' 

for any > 0. We compute in Fourier space, using that r/(p) = L^riQ{pL) (with an appropriate 
convention about the 27r), and f l6.69p . we get 

\\f{-zeV)vg{-tm\h= J I dpdq\f{ip)\'\fi{p-q)\'\g{iq)\' 

= {L/lf j I dpdq\f{p)\'\rfo{ip-q)L/i)\'\g{q)\' 
<CM{L/tf I [ dpdq 



<CN{l/Ld) 



2N 



L\6 



[l + {p-q)^{L/i)Y 
\l\\9\\l 



(6.73) 



which proves (16.701) . For the proof of fl6.7ip . we extract a decay of order [i/Ldy^ ^, estimate 
\g{(l)\ ^ llS'lloo and then integrate out q. The proof of f l6.72p is similar: 



wmvgmi'ns 
f{p) 2 



< C 



N 



1 +p2 



g{q){l + q 



"^KldJ 



dpdg 



:i+p') 



;i + g2)3[i + (p_g)2] 



and the last integral is finite. This completes the proof of Lemma [6.1 II 



(6.74) 



□ 
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6.2.4 Error terms in 06.421) for \i\ < Iq 

In Lemma [6. 101 we have estimated tlie terms i > iqiyi tlie last two summations in (16.421) . Note 
that the terms with i < —iq in the second line of fl6.42p identically vanish. It now remains to 
estimate the last four sums in fl6.42p for < |z| < zq- 

Lemma 6.12. Define Wi := 2^^^w and assume (16.261) and LW^^'^ <^ w. Then for any 5 > 
and h < hs we have 

F^H{A)7]-f7]+ J2 ^T^r F^H{A)7]{-f -l)r] 

0<i<io — «()<j<0 

^ E S^r [D^i^V^] + ^^T^r[D^FMl-l)v] 

0<i<io — io<j<0 

-C5h-^W^I^ f |V®A|2 (6.75) 

J B(2L) 

and 



Y F^H{A)riPri + ^ ^^Ti Ff H {A)r]{P - l)ri 

0<i<io — io<i<0 

<C6h-V/' [ \V^A\' + C^W'/'(—^Y. (6.76) 

Proof of Lemma \6.12[ We will present the proof for the summations over < i < io, 
the estimate of the negative i's are identical. We start with the first term in (16. 75 p . Since 
Wi = 2l*lw, then Wi < 2 for \i\ < io. Since on the support of Fj, > 0, it holds that 
^WiW <D^-W< \wiW, we obtain 

n 3 o 3 WjD^F? Wi , o , o 

> ^M^Ff > > ^^(D^ + W)Fi, 0<^< 

The analogous estimate for negative i will be 

HoF^ < -I^jWF^ < -^{D' + W)F^, 0>^> -i^. 

The additional W is necessary only for i < 0, when Tr Firj'yriFi may not be comparable with 
Tr D^FiTj'yriFi, but it is always comparable with Tr {D"^ + W)FiT]'yr)Fi. 
Using the second identity in (I6.38p . we have for < i < io, 

3fJTr F^H{A)T]-fr] =Tr F^HoT]-fr] + 23?Tr F^DA'r]-fr] + /i3?Tr F^a{B)r]-fr] 

ID ' 

>-^Tr (D^ + W)Fir]-fr]Fi - 2|Tr F^DA'r]-fr]\ - h\TT F^a{B')r]-fr]\. (6.77) 
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The last term was already treated in fl6.46p even after summation over i. 

To control the DA' error term, we proceed as before in fl6.48l) and fl6.6ip by inserting 
Gi:=l- h 

|Tr FfDA!r]-ir]\ < |Tr F^DA'FiT]-fT]Fi\ + |Tr F^A'Gir]-fr]FiD\. (6.78) 
The first term is estimated similarly to fl6.49p : 

1/2/ _ ~\l/2 



1/2 



|Tr F^DA'FifjjrjFil < (tt [AfFiTjjrjF,^ (tt FfO^r^^r^F, 

< c(Tt [AfF^y^^(TT D^F,7]-f7]F,y^^ 

< c(h~''^WiW''^/^ J [Afy^^ (Tt D^FiTi-iTiFi 
<SwiTT D^Fir]-fr]F, + C6'\h-^L^W)h-V/^ [ |V ® Ap. (6.79) 

Jb{2L) 

Here we used FiFi = Fi to remove the tildes from the terms with D^. We also used f l6.32p 
and the estimate 

snp Fi{x,x) < sup Fi{x,x) < Ch~^WiW^''^. (6.80) 

X X 

The second term in (16.780 is estimated similarly to fl6.65p by using Lemma 16.91 and the 
fact that ||V/i||i < CW. 

~ \ 1/2 / \ 1/2 



|Tr F^A'Gir]-fr]F,D\ < (Ti F^A'G^^A'F^^ (Tt D^F^t]-/!]'^ 

/ r \ 1/2 / \ 1/2 

<C[Lh-^W |V®An [Tt D^Fi7]-f7]Fi 

^ -'b{2L) ' ^ 



ru/1/2 r 

< SwiTr D^FiTi-friFi + C^"^ h'^W^^^ / |V ® Ap. (6.81 

Wi Jb{2L) 



Summing up fl6.79p and fl6.8ip for i < Iq < C\\ogh\ and using that LW^/"^ <^ /i^/^|log/i| ^ 
and LW^/"^ <^ w, we obtain 



^ |Tr F^DA'r]-fr]\ < 5 w^Tr D^Fir]'jr]Fi + Sh-V^^ j | V ® A| 

0<i<jo 0<i<io JB(2L) 



(6.82) 



if h < hg. Combining ( I6.82p with the positive terms from fl6.77p and choosing 6 sufficiently 
small, we obtain 

3? V Tr F^H{A)r]-fr] > V ^Tr {D^ + W)Firi-friFi - 5h-^W^/^ I |V ® A\^ . (6.83) 



0<i<io 0<i<jo 
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This completes the proof of (16.751) . The proof of fl6.76p is analogous by using that the mo- 
mentum supports of Fi and P are separated, so one can apply (I6.72p similarly to fl6.68p to 
obtain 

Tr FfHoTiPr] < Tr DFir]Pr]FiD = \\DF,7]P\\ls < CnW'/' { ^-^^y^ J ■ (6-84) 

The details will be left to the reader. This completes the proof of Lemma 16.121 □. 

Inserting the estimates from Lemma 16.81 Lemma 16.101 and Lemma 16.121 into fl6.42p . we 
have proved the following 

Proposition 6.13. Under the conditions of Theorem \6.4\ and assuming LW^^"^ <^ w, for any 
6 > and h < hs we have 

Tr HiA)T]^ri > Tr ri[HoU - C5^' [h + w^] W'/^Al - C^W"^/^ ( ^^^,^, ) 

+ [D'F^T^jT^] + J2 in^Tr [D'F^r^jr^] + ^ ^Tr [D^F^l " 7)^] 

0<j<io —io<i<0 

-C6h-^W^/^ [ |V®A|2. (6.85) 

Jb{2L) 

Finally, given < a < 1 and e > as in Theorem 16. 4| we choose w = h"^'^ and an integer 
> 3e~^, then Proposition 16.131 implies Theorem 16.41 □ 

As we already mentioned. Theorem 16.41 implies Theorem 16.11 

Upon combining Proposition 16.131 with the upper bound we get bounds on the kinetic 
energy terms in (I6.85p . These bounds will be used to control the substitution of Aq by Aj. in 
Section 16.31 

For each fixed u E B{^Lq), we will use Theorem [6]4] with the constant potential W replaced 
with Vu defined in (I6.1ip . We know that Co/2 < Vu < C , so factors W = can be replaced 
by constants in the estimates. To indicate the w-dependence, we define := — Vu- Recall 
that the operators Fi and F> defined in Section [6.2.11 depend on W . We will denote them by 
Fi and F" in case of 1^ = K- 

Theorem 6.14. Recall that Lq = and a = 1 — 3eo with Eq sufficiently small. Choose 

S{A) < £{0) + CH^-'^^'^Alo, (6.86) 

then, 

(111 

r.-3 <C/ii"^^°Aio, (6.87) 
-^1 
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where Tu is given by 



0<i<jo 



-io<i<0 



where 



lu ■= l(-oo,0] {Xu4>[Th{A - Cu) - Vu](t)Xu)- 

with Cu defined in (\6.9\\ and we set rju := Xu4'- 
Proof. By fl6:T0ll 



(6.88) 



(6.89) 



Tr 



<P[T,{A) - V]4> 



> / Tr 



Xu(p[Th{A - Cu) - K]0X« 



du 



(6.90) 



with Vu from (16.111) . Since V{0) = Cq by (15. 9p and using Lq <C 1 and the bound on W, we can 
assume that V > cq/2 on B{jLq). Also, using (I6.86P we get as in the proof of Corollary 16.31 
that 



1 



2Kh^ 



IV® < ch 



A 



B{2Lo) 



Therefore, (using the smallness of Eq) the condition fl6.14p is satisfied for all u in the domain of 
integration. Note that this is exactly the condition (I6.28P in Theorem 16.41 with W = Vu, with 
L = Li and shifting the center of the ball B{2L) to u. The condition (I6.27P is satisfied since 
we consider the vector potential A — Cu- Finally, the condition (16.260 is satisfied by (I6.13p . 
We can thus apply Proposition 16.131 and we therefore have 



Tr 



<p[n{A)-v]<p 



> 



Tr xu^o - -Cs[h + w^]V^^''Al, 

N 



h 



.1/2 



Estimating as around (I6.17P (and choosing sufficiently large), we get 



(6.91) 



Tr 



<P[n{A) -V]<P > I Tr Xu<P[Ho - V{u)]^<Pxu% - C6h~' [ | V ® A\' 

^1 JB(2Lo) 



- (C«r + Ch'^Lf + Cltj^L, + j r„^. (6.92) 



46 



We evaluate the main term using precise semiclassics. Define = 0(Lox), V{x) = V{Lqx). 
Then we get by unitary scahng and applying Theorem 15.61 with the effective semiclassical 
parameter h := H/Lq = h^/'^+^o \]xgx 

£(0) = Tr (^[-/i^A - = Tr {4>[-h^A - V]4>) _ 

\2\J2 



2{2'Kh)'^ / (j){xY[p' -V{x)]_dxdp+0{h' 



x) \P - V{x)]-dxdp + 0(AL„r). 



{2'Khf 



By Remark 16.21 



Tr Xu4>[Hq - 1/(M)]_(/)y„,^ = ^r^h~ 



V 



5/2,2 



0{Lih 



Inserting the choices of Lq = h^/"^'^^ and Li = /i^/^+^o gg^^ 



£{0)- j lTXumo-V{u)Y 



dn 



(6.93) 



(6.94) 



(6.95) 



Finally, choosing sufficiently small 5 (depending on and using (16.861) and (16.951) . we get 

dnizD. □ 



6.3 Estimates on the smoothed vector field 

This section is a technical preparation for the next Section Let A G if^(M'^) be a vector 
field. We fix a radial function x ^ C'^{B{1)) with Jx = 1 and define Xr{x) = 'r~^x(a;/r) and 
Ar = A* Xv We note that V ■ = on a ball B{L) if V ■ A = on B{L + r). 

Suppose also given cj)' with < 0' < 1, supp0' C B{3L/2), 0' = 1 on i?(L) and |V0'| < 
C/L. With this function (p' we define A' = (p'A and = then A = A' and = A'^ on 
B{L). The constants in the following sections may depend on x and on the constant C in the 
estimate \V4>'\ < C/L, but we will neglect these dependences. 

Lemma 6.15. If r < L/2, we have 

\\A'-AX<Cr^[ |V®A|2 (6.96) 

Jb{2L) 

and 

\\A'-A',\\l <C [ |V® (6.97) 

Jb(2L) 
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Proof. We have 



\A'-AX< 



\A-Ar\^ = 
B(3L/2) Jb{3L/2) 



B{1) 



[A{x) — A{x + rz)'jx{z) 



dx 



Air . 
< — 



B(3L/2) 



/ / rz ■ V Aj{x + trz)x{z) dtdz 
Jb{i) Jo 

y2 f f \z -VAjix + trz)\'^ dtdzdx 

J 3(1) Jo 



(6.98) 



by Cauchy-Schwarz. Upon changing variable in the x-integral, we obtain f l6.96p . 
By the Sobolev inequahty, we have 



A' - A'X <C / \V®{A'- A',)\' < C 



B(3L/2) 



|V ® (A - Ar)\^ + L~^\A - Ar\^ . (6.99) 



In the first term we use Young's inequahty, /g(3^/2) 1^ ® = Jb(3l/2) * (V ® A)p < 
Ib(2L) ® ^1^' ^fter extending the integration to MP and using r < L/2 to insert a cutoff 
function 1b{2L)- In the second term we use the previous calculation (16.981) . □ 



Lemma 6.16. Suppose the vector field A satisfies 



A = 0. 



(6.100) 



B{2L) 



Then, for all r < L/2 we have 



\Ar\\'L2(^B{3L/2)) ^ ^ 



IV® A| 



(6.101) 



B(2L) 



and for all multi-indices n G \ {0}. 



r\\L^{B{3L/2)) 



< C„r'-'l"l / |V 

J B{2L) 



A\ 



(6.102) 



where the constants C,Cn are independent of A,L and r. 
Proof. Analogously to the proof of Lemma 16.151 we get 



(6.103) 



B{2L) 
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Also, using f l6.100p and the Poincare inequality, 



\MmB(3Lm) < \\MmB(2L)) < CL^ / | V ® A\^ , (6.104) 



B{2L) 



where, in the first step, we used r < L/2 and Young's inequality as above. This finishes the 

proof of mm) . 

To prove the estimate for the derivatives, for the given multi-index n G N^\ {0}, we choose 
e G with lei = 1, such that n' := n — e & N^, i.e., Inl = In'l + 1. Then calculate 



L2(B(3L/2)) 



B(3L/2) 



{d-x){yWA){x-ry) dy 

B(l) 



2 



dx 



<r^-'l"l||5"'xll2 / / \{d^A){x-ry)\'dxdy 

Jb{1) Jb{3L/2) 

<Cr2-2|"l f \V®A\^dx, (6.105) 

Jb(2L) 

where we used the assumption that r < L/2, so B{3L/2) + B{r) C B(2L). □ 

We will use the results of this section for L = Lq = h^/'^^^° where Sq is a small positive 
number. First, with the choice of L = Lq, by change of variable and using r < Lq/2, we 
get Lemma 15.71 as a consequence of Lemma 16.161 Second, we will use as an input the apriori 
bound (15.1 op , which was already proven as a result of Corollary 16. 3[ namely 

/i-M |V® < C/i°Ai^, a = l-3eo- (6.106) 

Jb{2Lo) 

This will give bounds on various norms of A' — A'^ and d'^A^ that will be used in the next 
section. 



6.4 Smoothing A 

In this section we refine the result of Theorem 16.11 and complete the proof of Theorem 15.21 
by proving (15. lip . The vector potential A will not be removed as in Section [6l2| but rather 
replaced with Aj. and this results in a smaller error. We fix three length scales r, 
with 

Lo = h^/^-'\ Li = h^'^+'\ r = h^/^+P, (6.107) 

where p will be chosen as p = Ceq with a sufficiently large constant, and we assume that 
A satisfies (I6.106p . We will perform a similar analysis as in Section 16. 2[ in particular, we 
will again perform a dyadic decomposition in energy and the parameter w in the dyadic 
decomposition will be chosen as 

U7 = ^ /i5-i^o_ (6.108) 
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Lemma 6.17. For any sufficiently small p > there exists a positive constant Eq (in fact 
^0 = cp can be chosen where c is a universal positive constant) such that the following is 
satisfied. Let Lq and r be given by fl6.107p . and G C^{B{Lq)) with |9"</)| < C„Lq Let V 
satisfy the assumptions given in Theorem \5.2[ Let furthermore, A satisfy fl6.106p andV-A = 
in B{Lq + r). Then, 

Tr [0[T,(A) - V]cf\ _ > Tr [0[T,(A.) - V]ct>] _ - C/i^+^«Ai„. (6.109) 

Since f lS.lOp was already proven in Corollary I6.3[ Lemma 16.171 implies (15. lip and therefore 
finishes the proof of Theorem 15.21 The rest of this section will be devoted to the proof of 
Lemma 16.171 

Proof. We will use localizations as in Section 16.11 with L = Li. Notice that the apriori 
assumption (I6.14p will now be satisfied on all the small boxes Bu{2Li) with u G B{^Lq), since 

h-^ [ \V®A\^<h-^l |V® A|2 < C/i"(^)^Ai, = C/i1-9^«Al,, (6.110) 

Jbu{2Li) JB{2Lo) ^LiJ 

where we used (I6.106P and the fact that Vu > Cq/2. For the Pauli case, Th{A) = [a - (D + A)]"^, 
consider 

ThiA) = ThiAr) + aiA - A,)a(i?) + aiD)a{A - A) -{A- Ar) ■ (A + Ar) 

= Th{Ar) + 2D{A - Ar) - (A - Ar) ■ {A + Ar) + h(r{B - Br). (6.111) 

Here we used that A and Ar are divergence free which property holds on the support of 0. 
Note that (I6.11ip will be used only on supp (p. 
We denote 

7o:=l{-oo,o)(0[T,(A)-\/]0) (6.112) 
and we collect a few apriori estimates on •jq. 

Lemma 6.18. With the definition (I6.112P and assuming that A satisfies (16.1060 and V is 

bounded, we get 

Tr 7o < C/i"9^"/^Ai, (6.113) 
Tr ^^700 < CH-^'^/^Alo- (6.114) 
Proof. By the variational principle and Th{A) > {D + A)"^ — \B\, 

Tr 70 < 2Tr lHr:^)1^^oo,o) (0( (D + A)^ - \B\ - V)(t>) 

< 2Tr z.2(M3)l(_oo,o)((^ + A)^ - {\B\ + V)lsupp<^). (6.115) 

Notice that the last inequahty uses the fact that we consider the strictly negative eigenvalues. 
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By the CLR estimate we get, since V is bounded, 

Tr 70 < Ch~^ I {\B\ + \V\f'^ < CAlo + Ch~^ j \Bf^ 

Jsnmd) JB(Ln) 



suppt/) J B(Lo) 



< CAl, + Ch-'Ll'^K'^-KL.f = Ch-'^-l'k:^, 



using fl6.106p and a = \ — S^o, which proves fl6.113p . 
For the rest of this proof, set 

A' = A0', V' = V(I)\ (6.116) 

where </>' is a cutoff function such that supp</)' C B{3Lq/2), (p' = 1 on supp0 = B{Lq) and 
|V0'| < C/Lq. By the magnetic Lieb-Thirring inequahty. Theorem 12 .4^ we also have, 

>Tr (P{Th{A') - r)07o > "C/^"" / i^'f^ - / |V ® ^?)'^'( j [^f )'^' 

> - CAl, - Ch-^Lt{h''^'^KL}j"' = -Ch-'>'^^/'ALo, (6.117) 

using fl6.106p and a = 1 — Seq. Therefore, by Schwarz and Lieb-Thirring inequahties, from 
fleruTI) we get 

Tr ^20700 <2Tr [Th{A') + [Af ]07o0 < 2Tr [V + [Af ]07o0 + Ch~^'"/'ALo 
<2(| [V'+[Aff'f\l {M4>f'f\ch-'^o/^A,, 

<C^Ll + P'llf P'llf )'^'(/i-2Tr DW)'^' + C/i-9^"/'Alo, (6.118) 
where Qq is the density of 79. Since /b(2Lo) ^ — ^ (^^e fl5.7l) ). we have 



and 



\A'\\i< A'<CLi \V®A\' (6.119) 

'b{2Lo) Jb{2Lo) 



\A'\\l< [ \V (^A'\^ <C [ [\V^A\^ + Lo^\A\^]<C [ \V A\^ (6.120) 

J J B(2Ln) JB(2Ln) 



'B{2Lo) JB{2Lo) 

by Poincare and Sobolev inequahties, thus 

Jb(2Lo) 
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Clearly L^'^^^^h ^ L^^^ , so we can continue f l6.118p as 



from which (16.1141) follows since h ^^°/^Al(, ^ 1 if is sufficiently small. 



□ 



Returning to the decomposition (16.1 lip , we first show that the effect of the quadratic term 
(in A) is negligible: 

Lemma 6.19. 

T:t{A-Ar)iA + Ar)(f)-fo(t><Cir/Lo)^/^^{h-^Ll)h-^-^'"[ |V ® v4|'. (6.121) 



In particular, assuming (16.1061) and that p > IOOeq, we have 



TTiA-Ar){A + Ar)(P^O 



(6.122) 



Proof. Analogously to the proof of Lemma r6.18[ we use Holder and Lieb-Thirring inequalities. 
The prime on the A and Ar denote localizations to B{3Lo/2) as in (I6.116p . Then, using 
\p and (I6.119p - (I6.120I) . we have the following estimate: 



I^'IIp < 



(yl-^)(A + ^)0po0 <{J [iA'-A',)iA' + A'X''\ 

VlOii A/ A> ||9/10|| , ./ nl/lO 



5/212/5 



(0PO0) 



5/3 



3/5 



< \\A' - AX'^'WA' - KWe'^WA' + AX^'^'WA' + X\\T'[h-'Ti Z^W} 



3/5 



'B(2Lo) 



1/10 



IV® v4| 



B{2Lo) 



9/10 



^-6/5-27eo/20^3/5 



< C(r/Lo)'/'°(/i-'L2)/i-2-2^o / |V®Ap 

'B(2Lo) 



(6.123) 

which gives (16.1210 . Here we also used Lemma [6.151 with L = Lq and Lemma [6.18[ □ 



The main problem is to estimate the current term in (16.1111) . i.e. the term 



Tr 



a{A - Ar)a{D) + a{D)a{A - A,.) 



Tr 



2D{A- Ar) + ha{B ~ Br 



(6.124) 



where we used that A and Ar are divergence- free on supp (p. We will apply localizations as 
in Section 16.11 with L = Li on this term. Since this is a first order operator, we can do so 
without localization error, using 



Li 



2 '^^ 

L2./ ^"Zfj 
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for the partition of unity defined in fl6.ip . So we have 



TiDiA- A,) 070 



Tr 



D{A - Ar)7]ulQVu 



(6.125) 



where we set rju '■= Xu<P- Similar expression holds for the terms in the first line of f l6.124p 
Thus, we have to estimate (see f l6.116p for notation), 

1 du 



Tr 



(6.126) 



{2D{A - a;) + ha{B' - B',))r]uioVu 

or, equivalently, 

J Tr [{aiA' - Al)a{D) + a(D)a(A' - Al))rj^^oVu 
Here we have set B' := V x A', 5^ := V ® A'^- We recall that, similarly to f l6.45p . we have 

[Bf+ f <c [ |v® 

J J B{2Ln) 



1 du 
LI' 



(6.127) 



We define, with Vu from (16. lip . 



and 



P=l(_oo,0](i^") 



Furthermore, we let F^, F^, F^, F^ be defined as in Section [6. 2. II with W = Vu- The parameter 
w entering the definition of the F's has been fixed in f l6.108p above as w = = /i5~2^o. 

Lemma 6.20. Suppose p > lOOeo- Then, for any fixed u G i?(|Lo) and with rju '■= Xu4' '"'^ 
have 

Tr [{a{A - Ar)a{D) + a{D)a{A - Ar))vuloVu] > -CH'+'^Al, - Ch'^'^Tu^ (6.128) 

where Tu was defined in (16.880 . 

We postpone the proof of Lemma 16.201 and first finish the proof of Lemma 16.171 Using 
(EHH), fl6TT2D . the estimate (Kim from Lemma EH and fl6T25ll . we have 

Tr [0[T,(A) - V]<P]_ = Tr [<P[n{A) - V]<p-fo] 

> Tr [<P[n{Ar) - V]07o] - CH'^'^Alo 

Tr ^{a{A - Ar)a{D) + a{D)a{A - Ar))r]uloVu\ 
We now apply Lemma 16.201 and Theorem 16.141 to get 

Tr [0[T,,(A) - > Tr [0[T,(A) - V]<P^o] - Ch^^'°I^L, - Ch^'^ J Tu^ 

> Tr [0[T,(A,) - V]<P]_- CH'+^'^Al,. (6.129) 
This finishes the proof of Lemma 16.171 □ 
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Proof o f Lemma 1 6. 2(A From now on we fix m G B{Lq + Li) and drop the u indices and 
superscripts for simplicity, i.e, we set 7 = 7o and rj = 77^. We rewrite tlie current, for eacli 
fixed u, by using (16.391) . as 



Tr {a{A- Ar)cr{D) + a{D)(T{A- Ar))r]-fr] 



= Tr ^{a{A' - Al)a{D) + aiD)a{A' - K))ri{^ - P)ri 

= 2TT[F^DiA'-A'Ml-P)v] 
+ 2 5^ Tr [F^D{A' - A'Ml " 1)^] " 2 J] Tr [F^D{A' - A'^^P - 1)^] 




+ 2 J2 Tr[F^D{A'-A[.Uj-P)r^] 



0<i<io 



+ hJ2^^[F'^iB'-B'M^-P)rj] 



+ Tr ^Fl{(r{A' - A',)a{D) + a{D)a{A' - A',))r]{-f - P)r]^. (6.130) 



Note tliat we used formula (I6.126P for all terms with i < io, while we used (16.1271) for i > iq. 
Notice also that the sums over negative indices can be restricted to — ?o < < 0, since the 
Fj's vanish for larger values of Another important observation is that the left hand side 
of (I6.130p is real, so it suffices to estimate the real part of each term. The estimates will be 
very similar to the estimates of the terms (I6.48p . (16.611) and (I6.78P obtained during the apriori 
estimates in Section 16.21 but A will be replaced with A — A,, and we will capitalize on the 
the factor ^ gained from the smoothing in (16.960 compared with the usual Poincare 
inequality (I6.30p . Now we treat each term in (16.1301) separately. 

Step 1: The Fq term. We write 



Tr [F^D{A - A;)r/(7 - P)r,] 

= Tr [F^D{J^ - K)FoV{l " ^)^^o] + Tr [F^D{j^ - K)Gori{l - P)vPo]. (6.131) 



Recall that Fq is slightly larger than Fq, in particular FqFq = Fq and Fq + Gq = 1 (see Section 
16.2.11 for the precise definitions). The first term is estimated as in (16.491) . with = Fqij^Fq 
and its density qq as 




Tr D^uo 



1/2 



(6.132) 
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We also used that W = Vu is bounded. By the bounds on ||po||oo and Tr D'^ujq < CTr uq < 
CwAl^ from fl6.50p and f l6.96p . we have 



|Tr [F^D{A' - A;)For/(7 - P)r]Fo]\ < Ci^h-Wr'AL, 



B{2Lo) 



1/2 



(6.133) 



using fl6.106p and the choice of parameters f l6.107p . f l6.108p . 

For the other term in fl6.13ip . similarly to f l6.5ip and f l6.52p . we have 

|lV[F2Z}(A'-A:,)Gor/(7-P)r/Fo]| 

< CTr [D^FoV^FoY^^Tt [F^{A' - A';)Gl{A - A';)F^Y'\ (6.134) 

We apply the estimate Tr [D^Fot^^Fq] < CwKl^ as before together with the first inequality of 
f l6.54p from the application of Lemma 16.91 to get 

|Tr [F',D{A' - A'^)Gov{l - P)riFo]\ < (c«^AL,r/i-l V ® A||i.(^(2Lo))) 

<^^i+i{2p-2ieo)^^^_ (6.135) 

Here we used (I6.96P and that 



[\V®iA'-A[.)\'< [ \V®{A-Ar)\^ + CLf f \A-Ar. 

J J B{2Lq) Jb(3Li/2) 



<(l + C(r/Lo)2) / |V®Ap, 

Jb{2Lo) 

to collect the A-terms. In summary, we have proved 

Tr [F^{A' - A[.) ■ Dr/(7 - P)r/] < Ch^+'°AL, < C/i^+^^Ai^ 
assuming p > lOOeo- 

Step 2: i>io + l. 

Lemma 6.21. Assume that f l6.106p is satisfied and that p > lOOeo- Then we have 

Tr [F^D{A - Ar)7]iv] < 2/i^/^Tr [D^Fyr]-fr]Fy] + CH^/^Al,. 

Proof. Recalhng E.>.o+i = ^1 F>F> = F> and F> + G> = 1, we get 

Tr [FlD{A - Ar)7]-f7]] 

= Tr [F^D{A' - A[.)F^r]-fr]F^] + Tr [F^DiA' - A[.)G^r]-fr]F^]. 



(6.136) 



(6.137) 



(6.138) 



(6.139) 
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The first term we estimate as 



Tr [F^DiA' - A'^)F^7]-f7]F>] < Tr [F^D^F>7]-f7]F>]TT [{A' - A'^YFy7]-f7]F: 



1/2 



(6.140) 



In the last factor we apply Holder, Sobolev and Lieb-Thirring inequalities similarly to f l6.123p 

to get 



Tr [{A' - A'^fFyr^-ir,F^] < \\A' - A'^Wl^'W - AX'^i^h-'Ti [D''Fy7]-f7]F^ 

< r^l^h-^'^ [ I V ® (Tr [D'^FyTi-i'qFy^ '^'^ 

Jb{2Li) ^ 



\i ||9/5 



3/5 



(6.141) 



Inserting this into (16.1401) . using < 1, we find 



Tr [FlD{A' - A',)F>r]-fr]F^ 



< ^1/10^-3/5 



N 1/2 / _ _ n4/5 

|V®An (Tt [D^Fyr]-fr]F^] 



B(2Li) 



= r^/^Tr [D^F>r]-fr]F^] + T-^r^/^h'^f [ |V ® A, . 



(6.142) 



for any r > 0, and where we used (16.1061) to get the last estimate. 

We now consider the second term in (I6.139p . By a Cauchy-Schwarz inequality and by 
applying Lemma [6.91 (recall (I6.33P ) similarly as in (I6.134p - (l6.136p . we find 



Tr [FlD{A-A;)G^r,^r,F^] < (Tt [D^F>(r/7r/)^F>]Tr [F^{A' - A',)GUA' - A',)F?- 

< rTr [L)^F>r/7r/F>] + r'^h^^r / |V A\^ 

J B{2Li) 

< rTr [D2F>r/7r/F>] + r^^/it+^-^^oAig, 



1/2 



where we used (I6.106P to get the last estimate. 

Combining (I6.142p and (I6.143p . we get (16.1380 by choosing r = h^/^. 



(6.143) 



□ 



Step 3: |i| < io + 1- 
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Lemma 6.22. If p > IOOeq, 

i<0 i<0 

+ 2 J2 T^r[FiD{A'-A',U^-P)v] 

0<i<io 

+ hJ^^^T [fMb' - B'Mi - P)v]} 

i<io 

i<io 

- Ch^'° J2 2'''wTr [D^F^7]{1 - 7)7/] - (7/i^+"«Al,. (6.144) 

0>i>«o 

Proof. Consider < i < iQ. We write 

Tr [F^D{A' - A[)7]-fri] = Tr - A'^)Fi7]-f7]] + Tr - A^)G,r/7r/], (6.145) 

using + Gj = 1. We can estimate the first term, using FjFj = Fj, as, 

Tr [F^D{A' - K)F,r]^v]\ < {Tr [F,'D'F,7]^r]Fi]y/^{Ti [{A' - A'^^F^T^-fT^F,]}'/' 

< /i^"°2*wTr [D^F^r]-fr]] + h'^'''2-'w-^TT [{A' - A'^^F^] 

< h^'°2'wTr [D^F^ri-fT]] + Ch-^'°h-^r^ I \V ® A\^ 

Jb{2Lo) 

< h^'^TwTi [D^F^r]jr]] + C/i^+^^-^'^Alo, (6.146) 

where we used fl6.80p and (16.961) in the third inequahty. Then we used (I6.106P and the choice 
of parameters (I6.107P to finish the estimate. This is in agreement with the desired estimate if 
p > IOOeq (since the number of terms io in the sum is only logarithmic in h). 

The remaining term in (I6.145P is estimated using Lemma and that ||V/.j||i < CW < C: 

Tr [F,DF,{A' - K)G,'n^'n]\ < {Tr [F^D'r^^r]]^ {t^ [G^' - A',)F^{A' - A;)G,r/7r/] j'^' 

< h'^^ ■ TwTt [F^D'r^jr^] + V ® Ml^^^^,,^-^^ 

< /i^^o ■ 2*w;Tr [F^D^r]-fr]] + (7/i^+^"^"«Alo. (6.147) 

This is in agreement with (16.1440 if p > lOOeo- 

We also estimate the corresponding term with P. 

r ^ 1/2 

Tr [F^D{A' - A',)r]Pr]] < VTr [FfD''r]Pr]\Ti [Pr/(A' - A',)Ff{A' - A';)r,P\ 

1/2 



\DF,7]P\\hs[^^ [Pri{^ - K)Fh^ - K)VP] ■ (6.148 
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Using (I6.84P and that (h/wL) = h^°/'^^ the first factor can be made smaller than an arbitrarily 

1 /2 

large power of /i, while the second one is bounded by A^^ — y4J,||2, which is not bigger than 
a fixed positive power of h. So this term is negligible. The similar terms for negative indices 
i are estimated in the same manner with only notational changes. 

Also the a{B' — B'^) terms are readily controlled. We leave this part to the reader. □ 

We can now finish the proof of Lemma 16.201 We combine (16.1301) with f l6.137p and the 
results of Lemma 16.211 and I6.22[ This finishes the proof of Lemma 16.201 □ 
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